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ABSTRACT 

We study the origin and properties of "extra" or "excess" central light in the surface brightness profiles of 
cusp or power-law elliptical galaxies. Dissipational mergers give rise to two-component profiles: an outer 
profile established by violent relaxation acting on stars already present in the progenitor galaxies prior to the 
final stages of the merger, and an inner stellar population comprising the extra light, formed in a compact 
central starburst. By combining a large set of hydrodynamical simulations with data that span a broad range of 
profiles at various masses, we show that observed cusp ellipticals appear consistent with the predicted "extra 
light" structure, and we use our simulations to motivate a two-component description of the observations that 
allows us to examine how the properties and mass of this component scale with e.g. the mass, gas content, and 
other properties of the galaxies. We show how to robustly separate the physically meaningful extra light and 
outer, violently relaxed profile, and demonstrate that the observed cusps and "extra light" are reliable tracers 
of the degree of dissipation in the spheroid-forming merger. We show that the typical degree of dissipation is 
a strong function of stellar mass, roughly tracing the observed gas fractions of disks of the same mass over the 
redshift range z ^0-2. We demonstrate a correlation between the strength of this component and effective 
radius at fixed mass, in the sense that systems with more dissipation are more compact, sufficient to explain 
the discrepancy in the maximum phase-space and mass densities of ellipticals and their progenitor spirals. We 
show that the outer shape of the light profile in simulated and observed systems (when fit to properly account 
for the central light) does not depend on mass, with a mean outer Sersic index ^2.5. We also explore how 
this relates to e.g. the shapes, kinematic properties, and stellar population gradients of ellipticals. Extra light 
contributes to making remnants rounder and diskier, and imprints stellar population gradients. Simulations with 
the gas content needed to match observed surface brightness profiles reproduce the observed age, metallicity, 
and color gradients of cusp ellipticals, and we make predictions for how these can be used as tracers of the 
degree of dissipation in spheroid formation. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
galaxies: nuclei — galaxies: structure — cosmology: theory 



1. INTRODUCTION 



Thirty years ago, iToomrd (1 1977b proposed the "merger 
hypothesis," that major mergers between spirals could re- 
sult in elliptical galaxies, and the combin ation of d etailed 
observations of recent merger remnants ^ chweizer 1982; 
Lake & Dres"sie3 119861; iDovon et al.1 119941: iShier & Fischer . 



1998 



James et al.l 1 1 9991; iGenzel et al.l 1200 ll; iTacconi et al 



2003 basvra et al.1 120061 l2007t iRothberg & Josephl 12004 
2006ah an d e.g. faint shells a n d tida l features around 
elUpticals (Malin & Carte^ 11980. 1983; Schweizei| 119801; 
ISchweizeF & Seitzen 119921 ; ISchweizer ,1996) have lent con- 
sidera ble support to this picture (e.g. iBarnes & Hernquisl 
|1992|) . Furthermore, in the now established ACDM cosmol- 
ogy, structure grows hierarchically (e.g. White & Rees 1978|, 
making mergers an inescapable element in galaxy formation. 

However, it has long been recognized that purely dissipa- 
tionless (collisionless) mergers of stellar disks cannot explain 
the high mass and phase-space densities of nearby ellipti- 
cals (especially those of relatively low mass < L*), which 
are far more dense than local stellar disks o f the same mas s 
JOsriker 19 80; Carlberg 1986 ; Gunn_1987; Kormendvll989t) . 
Of course, spiral galaxies are not purely colUsionless systems, 
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but contain interstellar gas in addition to stars and dark matter. 
Furthermore, most ellipticals have old stellar populations, im- 
plying that their last gas-rich mergers occurred at z > 1, and 
therefore their progenitors may be high-redshift spirals. These 
were likely more dense than nearby spirals, and observational 
evidence (Erb et al. 2006) indicates they had even larger gas 
fractions than their present-day counterparts (/gas ~ 0.5, with 
some approaching /gas ^ 0.8-0.9). Because gas can radiate, 
it is not subject to Liouville's Theorem, and processes related 
to gas dynamics and star formation can reconcile th e high 
phase space densities of elli pticals relative t o spirals dGunnl 
11987; Lake 1989; Schweizer 1998). In detail, Hernquist et all 
([T99 3) estimated that ~ 10% of the stellar mass must be added 
in a compact dissipational component to account for the cen- 
tral densities of typical ellipticals. 

The possible importance of gas dynamics and triggered 
star formation in mergers is reinforced by obs ervations of 
ultraluminous infrared galaxies (ULIRGs) (e.g. iSoifer et al.l 
1984a b), which are always associated with mergers in the 
local Universe (iJoseph & Wrighll Il985t ISanders & Mirabell 
mH). The infrared emission from ULIRGs is thought to 
be powered by intense starbursts in the nuclei of these ob- 
jects, originatinghi_compactj_ce^^ of gas 
(e.g. IScovilleet al.1 [19861: ISargent et alJ [19871 [T989I). which 
will leave dense stellar r emnants (iKormendy & Sa ndersI 19921 ; 
iHibbard & YuiJ 119991: iRothbe rg & Joseph 2004). More- 
over, observations of merging systems and gas-rich rnerger 

D ovon et al.1 119941; 
as well as post- 



remnants (e.g.. Lake & Dressier 1 986, 
IShier & Fischerl Il998l: iJames et all Il99 
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starburst (E+A/K+A) galaxies ( lGotoll2005h . have shown that 
their kinematic and photometric properties are consistent with 
them eventually evolving into typical ~ L* elliptical galax- 
ies. The correlations obeyed by these mergers and remnants 
(e.g. JGenzel et a l. 200lVRothberg & Joseph 2006ab. and ref- 
erences a bove) are similar to e.g. the observed fundamental 
plane and iKormendvl (1 1 977bh relations for relaxed ellipticals, 
and consistent with evolution onto these relations as their stel- 
lar populatio ns age, as well as the c lustering and mass density 
of ellipticals (iHopkins et al.l |2007c). 

The link between these processes and the formation of 
ellipticals may be manife st in their surface b rightn e ss pro - 
files. Early work by e.g. 'Kormend^ (Il977al): iKind (1 19781) : 
lYoung etal.1 ([19 78): Lauer (1985) and IKormendvl (Il985ah 
(see IKormendvl Il987at for a review) showed that typi- 
cal elliptical surface brightness profiles were not as sim- 
ple as uniform r'^"^ laws. Typically, central profiles in- 
terior to ~ Ikpc deviate from r'/'* laws fitted to the en- 
velopes of ellipticals, falling both above and below the in- 
ward extrapolation of the r'/^ law. Contemporaneously 
with the discovery of black hole -host galaxy correlations 
(e.g.lEo rmendy & Richstone 1995; Ferrarese & Merrittl206ot 
iGebhardt et al.ii2000.) . Hubble Space Telescope observations 
of the centers of elliptical galaxies established that typi- 
cal < ellipticals exhibit central "cusps" - i.e. a con- 
tinued rise i n p ower-law lik e fashion towards small radii 
(lLauer etal.l[T99 1. 1992; Cr ane et all \l99^ iFerrarese et all 
1994 iKormendy et al.i il994t iLauer et al.1 \l99^ IKormendvl 
19991) . whereas the most massiv e ellipticals appea r to ex- 



hibit central flattening or "cores." IKormendvl (1 19991) demon- 
strated in a number of cases that these "cusps" appeared 
in some sense to be "extra" light, i.e. a distinct component 
above the inward extrapolation of the outer profile measured 
at large radii. With the co mbination of HST and ground- 
based photometry (iKormendy et al.l l2008l) . it now appears 
that this excess is ubiquitous in "cuspy" or "power-law" 
ellipticals, with mass ranges and spatial extents compara- 
ble to those expected from observations of ongoing merger- 
induced starbursts (Hibbard & Yun 1999; Rothberg & Joseph 
2004j) and nuraerical simulations (IMihos & Hernquistlll994al: 
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iFaber et aH ([1997) showed that the presence or absence of a 
cusp or power-law nuclear profile is strongly correlated with 
other, global properties of ellipticals - cusp ellipticals tend 
to be more rotationally supported, diskier, and have slightly 
higher ellipticities. They argued that these differences rein- 
force the idea that at least the cusp^ ellipticals are the direct 
product of gas-rich mergers, with dissipation forming the cen- 
tral "cusp" and giving rise to correlated kinematic and photo- 
metric properties. The central excess or "extra light" in these 
cases may therefore represent a distinct imprint of the degree 
of dissipation in the spheroid-forming merger 

Numerical modeling over the past twenty years has also 
indicated that gas physics and star formation play key roles 
in shaping elliptical galaxies and has elucidated the relation- 
ship of mergers to the various phenomena described above. 
The possible relevan ce of these addition a l proc esses was an- 
ticipated already by iToomre & Toomr j (Il972h . who asked 
whether mergers would not "... tend to bring deep into a 
galaxy a fairly sudden supply of fresh fuel in the form of inter- 

' There are some differences in the literature in the use of the term "cusp 
ellipticals." Unless otherwise stated, we will use it to refer to ellipticals with- 
out a central resolved core/flattening - i.e. "power-law" ellipticals. 



steUar material ..." iBarnes & Hernquis3 (1 19911 [19961) showed 
that tidal torques excited during major mergers excite rapid 
inflows of gas into the cent ers of gal axies, providing the fuel 
to power intense starbursts (iMihos & Hernquist 1994d, 199g) 
and to feed rapid bla ck hole growth (Di Ma tteo et al.l 120051: 
IHopkins et al.l l2005bl[ah . Gas consumption by the starburst 
and dispersal of residual gas by su pernova-driven wind s and 
feedback from black hole growth dSpringel et al.ll2005ah. cul- 
minating in a pressure-driven blast -wave (e.g. IHopkins et al.l 
2006tf,|HoEk ins & Hernquistl l2006'), terminate star formation 
so that the remnant quickly evolves from a blue to a red 
galaxy. Provided that the interaction involved a "major" 
merger,^ the remnant will resemble an elliptical galaxy, with 
the bulk of its mass on large scales made from pro genitor 
stars which experienced violent relaxation (e.g. B arnes|[r988l 
1992; Hernquist 1992, 1993b), and the dissipation-induced 
starburst appropriatel y boosting the concent r ation and central 



phase space density (iHernquist et al.l 1993|: [ Robertson et al.l 
2006b; Naab et al.1 [2006aL ICox et al.1 l2006bl)" Moreover, 



IMihos & Hernquistl (ll994a') predicted that this process should 
leave an observable signature in the surface brightness pro- 
files of remnants, in the form of an upwards departure from 
the outer 'd e Vaucouleuri (119481) r'/'^-law distribution in the 
inner regions: i.e. a central "extra light" above the inwards 
extrapolation of the outer profile. 

Understanding the processes responsible for establishing 
the structural properties of ellipticals and their correlations 
has the potential of revealing the formation histories of these 
objects. However, notwithstanding major observational and 
numerical advances, little effort has been made to use the ex- 
tra light content of ellipticals in this manner, and most studies 
have restricted their focus to determining whether or not some 
extra light component is evident. This owes largely to the ab- 
sence o f a detailed theoretical fram ework: while the original 
work of lMihos & Hernquistl(|19 94a^ predicted that such cusps 
should exist, a more refined treatment of star formation and 
feedback, along with better resolution, is required for more 
detailed interpretation and modeling. For example, owing to 
limited spatial and temporal resolution and a simplified treat- 
ment of star formation, the "extra light" profiles predicted by 
[Mihos & Hernquist (1994a) generally exhibited much more 
severe breaks than predicted by state-of-the-art simulations 
and seen in recent observations (see Appendix [Bji. There 
have been considerable improvements i n these areas in recent 
years ( e.g. Sprin gel & Hernqui st 2003": Sprin gel et al.ll2005bl: 
ICox et al. 2006c), and we take advantage of the se refin ements 
here, and in comp anion papers, fHopkins et all (l2008dl) (here- 
after Paper I) and IHopkins et al.l (l2008hl[al ). to study galaxy 
cusps or extra light in both simulations and observed systems. 
Our objective in this effort is to identify the existence and 
understand the origin of different components that contribute 
to the surface density profiles of ellipticals, their cosmolog- 
ical scalings and relevance for the formation history of such 
galaxies, and their implications for global galaxy properties. 

In this paper, we focus on the extra light in our simulations 

* In a major merger, tidal forces are sufficiently strong to drive nuclear 
inflows of gas and build realistic spheroids. The precise meaning of ma- 
jor merger in this context is bluiTed by a degeneracy between the progen- 
itor mass ratio and the orbit iHernquist 1989; Hernquist & Mihos 199^ 
IBo urnaud et al. 2005), but both numerical I Younger et al. 2008) and obser- 
vational (Dasyraetal. 2006; Woods et al. 2006) studies indicate that mas- 
sive inflows of gas and morphological transformation are typical for mass 
ratios only below ~ 3 : 1 . Unless otherwise noted, we generally take the term 
"mergers" to refer to major mergers. 
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and in known cuspy elliptical galaxies. In §|2]and §[3] we de- 
scribe our set of gas-rich merger simulations and the observa- 
tional data sets we consider, respectively. In § |4] we compare 
different approaches for fitting the surface density profile, and 
attempt to calibrate various methods in order to recover the 
physically distinct (dissipational versus dissipationless) com- 
ponents in merger remnants. In § |5] we compare our simu- 
lations with and apply our fitted galaxy decomposition to a 
wide range of observed systems. In §|6]we use these compar- 
isons to study how structural parameters of the outer stellar 
light and inner extra light component scale with galaxy prop- 
erties. In §|2]we examine how the existence and strength of 
the extra light component is related to galaxy structure, global 
shape and rotation, and show that it drives galaxies along the 
fundamental plane. We investigate how this extra light com- 
ponent influences and is related to stellar population gradients 
in ellipticals in §|8] Finally, in §|9]we discuss our results and 
outline future explorations of these correlations. 

Throughout, we adopt a JIm = 0.3, J^a = 0.7, Ho = 
70km s"' Mpc"' cosmology, and appropriately normalize all 
observations and models shown, but this has lit tle effect on 
our conclusions. We also adopt a lChabrierl (|2003|) initial mass 
function (IMF), and convert all stellar masses and mass-to- 
light ratios to this choice. The exact IMF systematically shifts 
the normalization of stellar masses herein, but does not sub- 
stantially change our comparisons. All magnitudes are in the 
Vega system, unless otherwise specified. 

2. THE SIMULATIONS 

Our merger simulations were performed with the paral- 
lel TreeSPH code GADGET-2 (Springel 20 05|), based on a 
fully conservative formulation ( Springel & Hernquisj l2002h 
of smoothed particle hydrodynamics (SPH), which con- 
serves energy and entropy simultaneously even when smooth- 
ing lengths evolve adaptively (see e.g., iHernquisl Il993al; 
lO'Shea et al.l l2005h . Our simulations account for radia- 
tiv e cooling, optional heating by a UV background (as 
in i Katz et all |1996i: iD ave et all 1 19991 although it is not 
important for the masses of interest here), and incorpo- 
rate a sub-resolution model of a multiphase interstellar 
medium ( ISM) to describe sta r formation and supernova 
feedback (Springel & HernquistI 12003). Feedback from su- 
pernovae is captured in this sub-resolution model through 
an effective equation of state for star-forming gas, en- 
abling us to stably evolve disks with arbitrary gas fractions 
(see, e.g. [S pringel et al. 2005li ISpringel & Hernquistll200l 
iRobertsonet al.li2006aiia) . This is described by the parameter 
^eos, which ranges from q'eos = for an isothermal gas with 
effective temperature of 10^ K, to q'eos = 1 for our full multi- 
phase model with an effective temperature ~ 10^ K. We also 
compare with a subset of simulations whi ch adopt the star for- 
mation and feedback prescriptions from iMihos & Hernquis3 
( fi994 a d. 1996), in which the ISM is treated as a single- 
phase isothermal medium and feedback energy is deposited 
in a purely kinetic radial impulse (for details, see, e.g. 
iMihos & HernquistI (1 1994cl) '). 

Although we find that they make little difference to the ex- 
tra light component, most of our simulations include super- 
massive black holes at the centers of both progenitor galaxies. 
The black holes are represented by "sink" particles that ac- 
crete gas at a rate M estimated from the local gas density and 
sound speed using an Eddington-limited prescription based on 
Bondi-Hoyle-Lyttleton accretion theory. The bolometric lu- 
minosity of the black hole is taken to be Lboi = e^Mc^, where 



= 0.1 is the radiative efficiency. We assume that a small 
fraction (typically « 5%) of Lboi couples dynamically to the 
surrounding gas, and that this feedback is injected into the 
gas as thermal energy, weighted by the SPH smoothing ker- 
nel. This f raction is a fre e parameter, which we determine as 
in Di Matteo et al.l (l2005l) by matching the observed Mbh - c 
relation. For now, we do not resolve the small-scale dynam- 
ics of the gas in the immediate vicinity of the black hole, but 
assume that the time-averaged accretion rate can be estimated 
from the gas properties on the scale of our spatial resolution 
(roughly « 20 pc, in the best cases). In any case, repeating 
our analysis for simulations with no black holes yields identi- 
cal conclusions. 

The proge n itor g alaxy models are described in 
ISpringel et al.l (l2005b), and we review their properties 
here. For each simulation, we generate two stable, isolated 
disk galax i es, ea ch with an extended dark matter halo with a 
iHernqu is"? 1*199( 1) profile, r notivated by cosmological simula- 
tions^^ avarro e t alJll99^ iBusha et"an l2005). an exponential 
disk of gas and stars, and (optionally) a bulge. The galaxies 
have total masses Mvir = Vvir/(10G//[z]) for an initial redshift 
z, with the baryonic disk having a mass fraction = 0.041, 
the bulge (when present) having nib = 0.0136, and the rest of 
the mass in dark matter The dark m atter halos are as s igned a 
concentration parameter scaled as in lRobertson et alj (l2006d) 
appropriately for th e galaxy mass and redshift following 
iBuUock et al.l (l200ll) . We have also varied the concentration 
in a subset of simulations, and find it has little effect on 
our conclusions because the central regions of the galaxy 
are baryon-dominated. The disk scale-length is computed 
based on an assumed spin parameter A = 0.033, chosen to be 
near the mode in the A distribution measured in simulations 
(IVitvitska et al.ll2002h . and the scale-length of the bulge is set 
to 0.2 times this. Modulo explicit variation in these parame- 
ters, these choices ensure that the initial disks are consistent 
with e.g. the observed baryonic Tully-Fisher relation and 
estima ted halo-galaxy mass scaU ng la ws (Bell & de Jongj 
.200 k iKormendv & Freemanll2004t iMandelbaum et al...2006l 
and references therein). 

Typically, each galaxy initially consists of 168000 dark 
matter halo particles, 8000 bulge particles (when present), 
40000 gas and 40000 stellar disk particles, and one black hole 
(BH) particle. We vary the numerical resolution, with many 
simulations using twice, and a subset up to 128 times, as many 
particles. We choose the initial seed mass of the black hole 
either in accord with the observed Mbh-c relation or to be 
sufficiently small that its presence will not have an immediate 
dynamical effect, but we have varied the seed mass to iden- 
tify any systematic dependencies. Given the particle num- 
bers employed, the dark matter, gas, and star particles are all 
of roughly equal mass, and central cusps in the dark matter 
and bulge are reasonably well resolved. The typical gravita- 
tional softening in our simulations is ~ 20-50pc in the < 
systems of particular interest here, with a somewhat higher 
^ 50- lOOpc in the most massive systems (yielding an effec- 
tively constant resolution 0.01 in terms of the effective 
radius). In Paper I and Appendix iBl we demonstrate that this 
is sufficient to properly resolve not only the mass fractions but 
also the spatial extent of the extra light components of interest 
here (although resolution may become an issue when attempt- 
ing to model the very smallest galaxies, with < lOOpc and 
L < 0.01 L*, as discussed in § O. The hydrodynamic gas 
smoothing length in the peak starburst phases of interest is 
always smaller than this gravitational softening. 
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We consider a series of sev eral hundred sim u lations of col- 
liding galaxies, d escribed in [Robertson et alJ (I2006bll3) and 
ICox et al.l(l20063bh . We vary the numerical resolution, the or- 
bit of the encounter (disk inclinations, pericenter separation), 
the masses and structural properties of the merging galax- 
ies, initial gas fractions, halo concentrations, the parameters 
describing star formation and feedback from supernovae and 
black hole growth, and initial black hole masses. 

The progenitor galaxies have virial velocities V^i^ = 
55,80,113,160,226,320, and 500kms-i, and redshifts z = 
0,2, 3, and 6, and our simulations span a range in final 
spheroid stellar mass M, ^ 10^- lO'^M©, covering essen- 
tially the entire range of the observations we consider at all 
redshifts, and allowing us to identify any systematic depen- 
dencies in our models. We consider initial disk gas fractions 
by mass of /gas = 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 (de- 
fined as the fraction of disk baryonic mass which is gas) for 
several choices of virial velocities, redshifts, and ISM equa- 
tions of state. The results described in this paper are based 
primarily on simulations of equal-mass mergers; however, by 
examining a small set of simulations of unequal mass merg- 
ers, we find that the behavior does not change dramatically for 
mass ratios to about 3:1 or 4:1. The mass ratios we study are 
appropriate for the observations of ellipticals used in this pa- 
per, which are only formed in our simulations in major merger 
events. At higher mass ratios, the resul t is a small bulge in 
a still disk-dominate d galaxy (see e.g. lYounger et al.l l2008t 
iHopkins et al.ll200"83g l), which we do not study here. 

Each simulation is evolved until the merger is complete and 
the remnants are fully relaxed, typically ~ 1 -2Gyr after the 
final merger and coales cence of the BHs. We then analyze 
the remnants following ICox et akl (l2006bl) , in a manner de- 
signed to mirror the methods typically used by observers. For 
each remnant, we project the stars onto a plane as if observed 
from a particular direction, and consider 100 viewing angles 
to each remnant, which uniformly sample the unit sphere. 
Given the projected stellar mass distribution, we calculate the 
iso-density contours and fit ellipses to each (fitting major and 
minor axis radii and hence ellipticity at each iso-density con- 
tour), moving concentrically from r = until the entire stellar 
mass has been enclosed. This is design ed to mimic obser- 
vational isop hotal fitting algorithms (e.g. iBender et al.lll987t 
lBendeij[T988 ). The radial deviations of the iso-density con- 
tours from the fitted ellipses are expanded in a Fourier series 
in the standard fashion to determine the boxyness or diskyness 
of each contour (the 04 parameter). Throughout, we show pro- 
files and quote our results in terms of the ma jor axis radius. 
For further details, we refer to lCox et al.l (12006b) . 

We directly extract the effective radius as the projected 
half-mass stellar effective radius, and the velocity dispersion 
a as the average one-dimensional velocity dispersion within 
a circular aperture of radius R^. This differs from what is 
sometimes adopted in the literature, where 7?e is determined 
from the best-fitting Sersic profile, but because we are fitting 
Sersic profiles to the observed systems we usually quote both 
the true effective radius of the galaxy and effective radii of 
the fitted Sersic components. Throughout, the stellar mass 
refers to the total stellar mass of the galaxy, and the dynamical 
mass Mdyn refers to the traditional dynamical mass estimator 

a^R^ 

Miy„=k—^, (1) 

where we adopt A; = 3.8 (roughly what is expected for a 
iHernquisl d 19901) profile, and the choice that most accurately 



matches the true enclosed stellar plus dark matter mass within 
Re in our simulations; although this choice is irrelevant as 
long as we apply it uniformly to both observations and simu- 
lations). When we plot quantities such as R^, a, and Mdyn, 
we typically show just the median value for each simula- 
tion across all ^ 100 sightlines. The sightline-to-sightline 
variation in these quantities is typically smaller than the 
simulation-to-simulation scatter, but we explicitly note where 
it is large. 

3. THE DATA 

We compare our simulations to and test our predic- 
tions on an ensemble of observed surface brightness pro- 
files of ellipticals. Specifically, we consider three sam- 
ples of cusp or extra light ellipticals and a compilation 
of remnants of recent gas-rich m ergers. The firs t is th e 
V-band Virgo elliptical survey of iKormendy et al.l (l2008h . 
based on the complete sample of Virgo galaxies down to 
extremely faint systems Mg ^ -15 inp inggeli et al. (1985b 
(the same sample studied in Cote et al. 2006; iFerrarese et aP 
,2006,) . ,Kormendv et al. ( 2008) co mbine observa tions from a 
large numb er of source s (includinglBender et"an[T9 88, 200% 
ICaon et al."1990', '1994"; 'Davis et al."l985'; 'Jedrzeiewski et al] 
1987; Jedrzeiewski 1987; Kormendv et al. 2005; Lauer 19851; 
iLauer et alJil995L i2005i: iLiu et al.ll2005l; iPeletier et aljil99o[) 
and new photometry from McDonald Observatory, the HST 
archive, and the SDSS for each of their objects which (af- 
ter careful conversion to a single photometric standard) en- 
ables accurate surface brightness measurements over a wide 
dynamic range (with an estimated zero-point accuracy of 
±0.04 y mag arcsec"^). Typically, the galaxies in this sample 
have profiles spanning ^ 12-15 magnitudes in surface bright- 
ness, corresponding to a range of nearly four orders of mag- 
nitude in physical radii from ^ lOpc to ^ lOOkpc, permitting 
the best simultaneous constraints on the shapes of both the 
outer and inner profiles of any of the objects we study. The 
profiles include e.g. ellipticity, 04/0, and^-z colors as a func- 
tion of radius. Unfortunately, since this is restricted to Virgo 
ellipticals, the number of galaxies is limited, especially at the 
intermediate and high end of the mass function. 

We ther e fore a dd surface brightness p rofiles from 
iLaueret al] (l2007ah , further supplemented by iBender et alj 
(119881) . I Lauer et al] (l2007ab compile V-band measurements 
of a large number of nearby systems for which HST imag- 
ing of the galacti c nuclei is available. These include the 
ILauer et al J l2005h WFPC2 data-set, the Laine et al. {2003 
WFPC2 BCG sample (in which the objects are specifically 
selected as brigh test cluster galaxie s from Postman & Laued 
(1995)), and the iLauer et air (ll995h and Faberetal. (199^ 
WFPCl compilations (see also lOuillen et al.ii2000i ; iRest et al] 
2001; Ravindranath et al. 2001). Details of the treatment of 
th e profiles and conv ersion to a single standard are given 
in iLauer et al.l (12007 ai) . The sample includes ellipticals over 
a wide range of luminosities, down to Mb ^ -15, but is 
dominated by intermediate and giant ellipticals, with typi- 
cal magnitudes Mg ^ -18. This therefore greatly extends 
our sampling of the intermediate and high-mass end of the 
mass function, but at the cost of some dynamic range in the 
data. The HST images alone, while providing information 
on the central regions, typically extend to only ~ 1 kpc outer 
radii, w hich is insufficient to fit the outer profile. ILauer et al] 
(2007a) therefore combine these data with ground-based mea- 
sur ements from a number of sources (see the references for 
the iKormendv et all (|2008[) sample) to construct profiles that 
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typically span physical radii from ^ lOpc to ^ 10-20kpc. 
Althoug h the composite profiles were used in iLauer et aP 
(l2007al) to estimate effective radii, they were not actually 
shown in the paper It should also be noted that there is no sin- 
gle criterion that characterizes galaxies included in this sam- 
ple, but they generally comprise luminous nearby ellipticals 
and SO galaxies for which detailed imaging is available. We 
emphasize that issues of completeness and e.g. environment 
are not important for any of our conclusions. 

We occasionally supplement the profi les from Lauer et alj 
(p007a) with addi tional profiles used in iBender et al. ( 19881 
119921 1 1993L[T994h . and in some cases subsequently updated. 
These are more limited: typically the profiles cover ~ 7 mag- 
nitudes in surface brightness, extending from ~ 30-50pc out 
to ~ a few kpc (typically ~ 3 kpc in low-luminosity systems, 
and ^ 15 kpc in the brightest systems, sufficient for accept- 
able, but not strong constraints on the outer profile shapes). 
However, the measurements are usually in each of the V, R, 
and / bands, and hence allow us to construct multicolor sur- 
face brightness, ellipticity, and 04/0 profiles. We use this to 
estimate e.g. the sensitivity of the fitted parameters and galaxy 
profiles on the observed waveband and on the quality and dy- 
namic range of the photometry. 

In various places, we compare our results from the study of 
these ellipticals to our results in Paper I from a study of local 
remnants of gas-rich merger rem nants ( Rothberg & Joseph 
12004 . For these objects, Rothb erg & JosephI (12004) com- 
pile /T-band imaging, surface brightness, ellipticity, and 04/0 
profiles, where the profiles typically range from ^ lOOpc to 
^ 10-20 kpc. These span a moderate range in luminosity (in- 
cluding objects from ^ -20 to Mk ~ -27, but with most 
from Mk ^ -24 to Mk ^ -26) and a wide range in merger 
stage, from ULIRGs and (a f ew) unrelaxed systems to shell 
ellipticals. As demonstrated in lRothberg & JosephI ( l2004l) and 
argued in Paper I, these systems will almost all become (or al- 
ready are, depending on the classification scheme used) typ- 
ical ^ ellipticals, with appropriate phase space densities, 
surface brightness profiles, fundamental plane relations, kine- 
matics, and other properties. For a detailed discussion of the 
modeling of these systems and the profiles themselves, we re- 
fer to Paper I (all of the results shown for these systems are 
derived therein). We show the results from Paper I here in 
order to test the continuity of merger remnant and (cusp) el- 
liptical populations. 

Because we are here specifically interested in extra light or 
cusps in observed ellipticals, and because the generally ac- 
cepted belief is that core ellipticals are not directly formed 
in gas-rich major merg ers but are subsequently modified by 
dry re-mergers ( see e.g. lFaber et alJI 19971: Ivan Dokkuml2005l: 
iBell et al] l2006h . we restrict our attention only to those el- 
lipticals which are confirmed via HST observations as being 
cusp ellipticals. We include all the confirmed gas-rich merger 
remnants, but note there are a small number of extreme unre- 
laxed cases for which sharp features in the surface brightness 
profiles prevented derivation of meaningful quantities (note, 
however, as shown in Paper I, that almost all of the objects in 
this sample are sufficiently well-relaxed at the radii of inter- 
est for our fitting). We exclude dwarf spheroidals, as they are 
not believed to form in rnajor mergers as are ell ipticals (e.g. 
lKormendvlfT985R Il987bt iKormendv et al.ll2008l) . and in any 
case they dominate at extremely low masses where our sim- 
ulations do not sample the population (they also predominate 
as satellite galaxies, whose effects we do not model). 

We also exclude SO galaxies (adopting the morphological 



classifications from iKormendv et al.l ( l2008l) and iLauer et al] 
(l2007al) . although it makes little difference exactly which 
classifications we consider). This is not because of a phys- 
ical distinction: observations suggest that these likely form 
a continuous family with the low-luminosity cusp ellipticals, 
and in fact a number of our simulated gas-rich merger rem- 
nants would, from certain viewing angles, be classified as SOs. 
However, in order to derive e.g. the parameters of the outer, 
violently relaxed profile and central extra light, it would be 
necessary to remove the contribution of the large-scale disk 
from the surface brightness profiles of these objects. Our two- 
component (outer dissipationless and inner dissipational) Ser- 
sic models (described in § HJi then become three-component 
fits, and the degeneracies involved with three independent 
components, even with our best data and simulations, are so 
large as to render the results meaningless. We have, however, 
re-visited all of the SOs in these samples in light of our re- 
sults, and find that they are, in all cases, consistent with our 
predicted and observed trends. However, it is too difficult to 
infer these trends directly from the SOs themselves without 
ideal disk subtraction. 

This yields a final sample of « 80 unique elliptical galax- 
ies, and « 50 confirmed remnants of gas-rich mergers. Most 
of the sample spans a range of three orders of magnitude in 
stellar mass, from < O.IM* to ~ lOM*, and a wide range in 
extra light properties. There is, of course, some overlap in the 
samples that define our compilation; we have ^ 300 surface 
brightness profiles for our collection of « 80 unique ellipti- 
cals, including (for many objects) repeated measurements in 
multiple bands and with various instruments. This turns out 
to be quite useful, as it provides a means to quantify error es- 
timates in fits to these profiles. The variations between fit pa- 
rameters derived from observations in different bands or made 
using different instruments are usually much larger than the 
formal statistical errors in the fits to a single profile. There are 
no obvious systematic effects (i.e. systematic changes in pro- 
file fits from y to / band), but as demonstrated in Paper I the 
effects of using different bands or changing dynamic range 
(from different instruments) can be complex, depending on 
the structure and degree of relaxation of the outer regions of 
a system. On the other hand, there are well-relaxed objects 
for which almost no significant change in the fits occurs from 
band to band. It is therefore useful to have multiple obser- 
vations of the same system, as it allows us to get some idea 
of how sensitive our fits are to differences in e.g. the choice 
of observed wavelength or dynamic range from instrument to 
instrument. 

In Table [T] we list the names and properties of our sam- 
ple ellipticals, including the relevant sources of photometry, 
stellar masses, absolute magnitudes, stellar velocity disper- 
sions and effective radii, ellipticities, isophotal shapes, and 
rotation properties. We have converted all the observations 
to physical units given our adopted cosmology, and compile 
global parameters (where not available in the original papers) 
including e.g. kinematic properties, luminosities, and black 
hole masses from the literature. We determine stellar masses 
ourselves in a uniform manner for all the objects, based on 
their total /T-band lumin osities and {B - V ) color-dependent 
mass-to-light ratios from lBell et alj (l2003l) . corrected for our 
adopted IMF. We have also repeated our analysis using stellar 
masses derived from a mean M/L as a function of luminos- 
ity or from fitting the integrated UBVRIJHK photometry of 
each object to a single st ellar population with the models of 
iBruzual & Charloi (l2003l) . and find this makes no difference 
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to our conclusions. 

Throughout, we will usually refer interchangeably to the 
observed surface brightness profiles in the given bands and 
the surface stellar mass density profile. Of course, stellar 
light is not exactly the same as stellar mass, but in Paper I 
and §|8] herein, we consider the differences between the stel- 
lar light and the stellar mass density profiles as a function of 
time, wavelength, and properties of the merger remnant, and 
show that the V and /T-band results introduce little bias (i.e. 
are good tracers of the stellar mass); the Sersic indices and 
extra light fractions fitted to the /T-band profiles of the simula- 
tions are good proxies for the Sersic index of the stellar mass 
profile and extra mass/starburst mass fraction, even close in 
time to the peak episode of star formation. 

Although we are not concerned about the absolute normal- 
ization of the profile (i.e. mean M/L), since we derive to- 
tal stellar masses separately from the integrated photometry, 
we must account for systematics that might be induced by a 
change in M/L as a function of radius. The results from our 
merger remnant sample (observed in /T-band) are, on average, 
robust in this sense, but they should be treated with care, es- 
pecially in the most extreme cases (namely the few LIRGs 
and ULIRGs in the sample), where younger stellar popula- 
tions may decrease M/L towards their center (see § |8]i. We 
emphasize though that many of these systems are much older 
and more relaxed (e.g. ellipticals with faint shells or tidal de- 
bris). The profiles in optical bands such as V require more 
care - when the system is very young (< 1 -2Gyr after the 
major merger-induced peak of star formation), there can be 
considerable bias or uncertainty owing to stellar population 
gradients and dust. However, once the system is relaxed, the 
optical bands also become good proxies for the stellar mass 
distribution. 

In fact, in Paper I and §[8] we demonstrate that once the sys- 
tem reaches intermediate age, the bias in e.g. B or V band is 
often less than that in K band, because systems tend to be both 
younger and more metal rich in their centers. In /T-band, these 
both increase L/M, leading to a (small) systematic bias. In op- 
tical bands, however, the two have opposite effects (younger 
age increases L/M, but higher metallicity decreases L/M), 
and they tend to mostly cancel. Since essentially all of our el- 
lipticals are older than this stellar population age (even in their 
centers), and they have been carefully vetted and either cor- 
rected for the effects of e.g. dust lanes in the sources or (where 
correction w as too difficult) excluded from our samples (see 
iKormendv' ;t al. 2008; Laueretal. 2007a), we are not con- 
cerned that significant bias might persist. Furthermore, com- 
parison of systems observed in different bands demonstrates 
that our conclusions are unchanged (modulo small systematic 
offsets) regardless of the observed bands in which we ana- 
lyze these systems. As has been noted in other works, most 
of these objects have weak color gradients, indicating little 
variation in M/L with radius. 

4. RECOVERING THE PHYSICALLY APPROPRIATE "EXTRA 
LIGHT" 

We would like to use the surface brightness profiles of 
merger remnants to estimate the contribution from extra light 
and, in particular, to infer the fraction of stellar mass that was 
formed in a compact central starburst. However, as noted in 
Paper I, the light profiles in our simulation remnants are quite 
smooth, even where the extra light fraction is large. This 
makes recovering the extra light component a non-trivial pro- 
cedure, which can be sensitive to the assumptions made in 



fitting (see also iNaab & Truiiiloll2006h . We discuss various 
procedures and their consequences in Paper I, but briefly re- 
view them to highlight the most important decomposition we 
will adopt. 

Figure [T] shows the surface density profile of a typical 
merger remnant from our simulation library, with a gas frac- 
tion of ^ 10%, which happens to provide a good match to 
several observed ellipticals. We begin by reducing the pro- 
file to the two most physically relevant components: the "pre- 
starburst" or "disk" stars, i.e. those formed in the rotation- 
ally supported disks before the final coalescence of the galax- 
ies, and the "starburst" stars, produced in the final, compact 
starburst. We operationally define this as stars formed within 
±125 Myr of the peak in the starburst star formation, but since 
the starburst is usually very distinct, changing this definition 
within reason makes almost no difference. We combine both 
the stars present in the stellar disks at the beginning of our cal- 
culations and the stars formed in the disks over the course of 
the simulations (the disk stars and pre-starburst stars in Fig- 
ure[T]of Paper I) into a single pre-starburst stellar population, 
because there is no robust physical distinction between the 
two. Not only do both populations experience violent relax- 
ation in the final merger/coalescence, but it is also arbitrary 
where (or at what time) we initialize our simulations, relative 
to the final merger, and therefore what fraction of the stellar 
mass forms before our simulation begins. 

We also neglect embedded stellar disks which can be 
formed by gas remaining after the merger; as demonstrated in 
Paper I, these contribute negligibly to the surface mass den- 
sity profile in even the most gas-rich merger remnants (and 
we are primarily concerned with the surface mass density pro- 
files, not the early-time optical/UV light profiles which more 
strongly reflect the light from new stars). This does not mean 
such disks are unimportant - indeed they are ubiquitous in 
gas-rich merger remnant simulations and contribute critically 
to the kinematics (especially the rotation and 04, see Cox et al] 
2006b), and observations suggest that low-level disks may be 
present in nearly a U cusp galaxies (perhaps all, give n pro- 
jection effects; see iFerrarese et al.lll994t iLauer et al.l 12005). 
However, especially given our exclusion of SO galaxies, these 
are not a significant component of the surface brightness pro- 
files. 

We showed in Paper I that the total profile of the system 
can be robustly represented as a two-component sum of two 
Sersic distributions. If we fit each known physical component 
separately, we find that the pre-starburst component follows a 
nearly exact Sersic law (here, n., = 2.85), while the inner com- 
ponent (given its shape and origin in a gas-rich, dissipational 
event) can be well fit by a lower-Wj ^ 1 law.^ Based on this 
behavior, we therefore adopt a two-component decomposition 
of the observed quantity, the total surface brightness profile. 
This is defined as a Sersic plus cusp or extra light model, with 
an outer component for the pre-starburst stars with a free Ser- 
sic index, and an inner component reflecting the starburst stars 
(with fixed = 1 when the shape of the starburst is not well- 
resolved, or free - albeit still generally low - when reso- 

^ Formally, we consider our simulation profiles outside of some multiple 
~ 3 — 5 times the resolution limit, or with a seeing con'ection appropriate for 
the comparison observed samples. We equally sample the profile in log r over 
a dynamic range extending to the largest radii in the observed samples (§[3). 
and weight each point equally assuming an intrinsic 0. 1 mag point-to-point 
variance in the SB profile (the typical magnitude of residuals fitting arbitrary 
splines to the data). We have varied these choices and find that our fits and 
conclusions are not sensitive to them. 
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Fig. 1 . — Upper Left: Surface mass density of a typical merger remnant from our simulation library (black), decomposed into stai's formed prior to the final 
merger (which are then violently relaxed; red) and stars formed in the dissipational starburst (blue). The Sersic index fitted to the pre-starburst component alone 
is shown, with the stellar mass fraction of the starburst component. Upper Right: Two-component (Sersic plus cusp or extra light) fit (inner exponential and outer 
Sersic) to the total light profile, with the Sersic index of the outer component and mass fraction of the inner component, and rms scatter (crp, in magarcsec"^) 
about the fit. Lower Left: Single Sersic function fit to the profile. Lower Right: Core-Sersic function fit. Our two-component, cusp plus Sersic function fit (top 
right) accurately recovers the profile of the violently relaxed component and mass fraction of the starburst component. The other fits give less intuitive results in 
this case. 

important, this fit, despite having no direct information about 
the physical components into which we decompose the bright- 
ness profile, recovers almost exactly the appropriate parame- 
ters for both components. The best-fit Sersic index («, = 2.83 
compared to 2.85) and effective radius of the outer or bulge 
component reproduce well those from fitting directly to the 
pre-starburst stellar population. Likewise, the inner or extra 
light component is a good match to the starburst component, 
and the fit recovers the extra light fraction accurately (5.1% 
compared to 6.1%; a smaller difference than reasonable un- 
certainties in our physical definition of the starburst compo- 
nent). Again, we emphasize that we have simply fit a function 
of the form in Equation ^ to the total surface brightness pro- 
file, ignoring our knowledge of the genuine physical break- 
down; but we find that we recover an accurate reflection of 
that true decomposition. 

Figure |2] shows the results of repeating this procedure for 
several hundred simulations; we directly compare the fitted 
extra light mass fraction and size (7?extra) to the mass frac- 
tion and size of the known physical starburst component, and 
find that the fitted components recover the physical values 
in the mean with a factor ^ 2 scatter This result is ro- 
bust with respect to e.g. the mass, orbital parameters, mass 
ratios, initial gas content, treatment of feedback and model 
for the ISM equation of state, and redshift of our simula- 
tions. We similarly compare the effective radius of the entire 
galaxy determined from these fits to that known from the di- 
rect integration of the profile. The fits recover R,, to within 



lution permits). We have studied this decomposition in Paper 
I, and show that it provides a good description of both simu- 
lated and observed merger remnant profiles, accurately sepa- 
rating the central light of younger stellar populations (as ob- 
servationally determined) from the light of older stars, which 
form a more extended distribution. The total surface bright- 
ness profile is then 

^ ^ "extra ^ J 

+/oexp|-/7„(--^— )'^"'|, (2) 

where .Rextra and Roatei are the effective radii of the inner 
(n' ~ 1) and outer (free n,) components (which we identify 
with the starburst and old bulge or pre-starburst components, 
respectively), /' and Iq are the corresponding normalizations, 
n' is the Sersic index of the inner (extra light) component 
(fixed n' = 1 where resolution limits apply) and is the Ser- 
sic index of the outer bulge or pre-starburst component. The 
constant b„ is the appropriate function of tis such that /?extra 
and /?oiiter correspond to to the projected half-mass radii. 

The upper right panel of Figure [T] shows an example of 
the outcome of such a fit. The resulting model of the sur- 
face density profile fits the simulation well, with a rms de- 
viation (cTp = (A/i^)'/^) of only cr^ ^ 0.15 mag (assuming 
ji oc -2.5 log/tot). This is comparable to the point-to-point 
variance in the profile of this simulation if we fit an arbitrary 
spline to the data, and thus reflects a genuinely good fit. More 



Hopkins et al. 



8 




0.1 1.0 0.1 1.0 

Wbura R^(starburst) [kpc] 




1 10 1 2 3 4 5 6 7 

Re (true) [kpc] ns(pre-starburst stars) 



Fig. 2. — Success of our proposed two-component empirical decompo- 
sition at recovering the known pliysical parameters of the galaxy starburst 
and pre-starburst (violently relaxed) components. Top Left: Mass fraction in 
the fitted "extra light" component /extra versus the known mass fraction of 
the physical starburst fth- Each point is the average across ~ 100 sightlines 
to a given simulation, although the sightHne-to-sightline variance is moder- 
ate (Si 0. 15dex). Different colors and symbols denote different initial disk 
gas fractions, orbital parameters, and merger redshifts (see key in Figure [741 . 
The fitted /extra recovers the physical /sj, on average, with a factor ~ 2 scatter 
(plotted error bar), and without any significant bias from any varied simu- 
lation parameters. Top Right: Fitted effective radius of the extra light com- 
ponent versus the projected half-mass radius of the known starburst stars. 
Again, the true values are recovered with a factor 2 scatter, independent of 
simulation choices. Bottom Left: Fitted effective radius of the entire galaxy 
versus Re from direct profile integration. The fits recover Re to better than 
'^0.1 dex. Bottom Right: Fitted Sersic index of the outer (dissipationless) 
component versus that fitted directly to the known pre-starburst component 
profile. The outer profile shape Jij is recovered to within ~ 0. 1 dex (the only 
significant bias is when there is a large disk in the remnant; but those cases 
are not relevant for this paper). The empirical decompositions do well at re- 
covering the known parameters in the simulations, across the entire range of 
simulation parameter space we have surveyed. 

^ 20-25%, given a dynamic range comparable to the obser- 
vations with which we compare. We also compare the Ser- 
sic index of the fitted outer component to that fitted directly 
to the known physical pre-starburst component and find the 
two agree within ^ 30-40%, without any systematic depen- 
dences except where the remnants have large embedded disks 
(biasing the fitted decompositions to lower n,). This gives us 
confidence that such an empirical approach can be used, in a 
statistical sense, to recover physically meaningful parameters 
describing the galaxy components. 

We note that, so long as we are fitting radii > 50 pc, we 
obtain similar results for a free inner component shape pa- 
rameter or fixed inner n'^ =1. In Paper I we demon- 
strate that the choice of «j = 1 for the inner component re- 
covers, on average, the correct physical mass fraction which 
participated in the starburst, and mitigates against degen- 
eracies in fitting to the outer Sersic profile (so that we re- 
cover the same value if we fit the two stellar populations 
independently). We therefore adopt this choice for our de- 
compositions in situations where we cannot reliably resolve 
the innermost extra light shape/stru cture (namely our simula- 
tion s and the obs erved samples of Roth berg & JosephI (120041) 
and iBender et a l. (1988)). We emphasize, however, that this 
choice for the inner component does not imply that this re- 
flects the true shape of the central extra light, which can 
be complex (see e.g. the range of central profile shape s in 
iLauer et al.ll200llCote et aIll2006tlKormendv et al.ll2008l) . In 
fact, caution should be taken when considering the central 
^ 30-50pc, which our simulations do not generally resolve 
(for most of the galaxy observations considered here, this cor- 
responds to ^ 0.5", a factor ~ 10 larger than the HST diffrac- 



tion limit). We include a detailed discussion of profile shapes 
as r ^ in Appendix IB1 

For now, we simply note that the quantities of interest 
here are well-converged in resolution studies, and our nu- 
merical tests (smoothing extremely high-resolution simula- 
tions over various seeing) find that n, = 1 choice is robust 
(in the sense that the mean properties are recovered similarly 
in either case). However, where the information is available 
in the central regions of th e galaxy (namely th e samples of 
iKormendv etall (I2008h and lLauer etall (l2007al) ). we find the 
best results using all of the information and freeing the Sersic 
index of the inner light component, to accommodate the real, 
observed shapes and structure in the inner light component. 
Despite these caveats, our treatment describes the starburst 
mass profile well where it is important to the overall surface 
density (from ~ lOOpc to ^ 1 kpc) and accurately recovers the 
total mass in the starburst component and its effective radius. 

We apply this formalism to the observed systems in §|5]be- 
low'^, and we present the results of our fits to each elliptical in 
our sample in Table [T] For sources with multiple independent 
observations, we define error bars for each fit parameter rep- 
resenting the ^ 1 (T range in parameters derived from various 
observations, typically from three different surface brightness 
profiles but in some cases from as many as 5 - 6 sources 
(where there are just 2 sources, the "error" is simply the range 
between the two fits). In many cases the different observations 
are comparable; in some there are clearly measurements with 
larger dynamic range and better resolution: the errors derived 
in this manner should in such cases be thought of as the typ- 
ical uncertainties introduced by lower dynamic range or less 
accurate photometry. 

In terms of direct comparison with our simulations, the data 
often cover a dynamic range and have resolution compara- 
ble to our simulations, provided we do not heavily weight the 
very central (< 30pc) regions of HST nuclear profiles. Ex- 
perimenting with different smoothings and imposed dynamic 
range limits, we find it is unlikely that resolution or seeing 
differences will substantially bias our comparisons. They do 
introduce larger scatter; the robustness of our results increases 
considerably as the dynamic range of the observed profiles is 
increased. 

As we demonstrate in Paper I, care should be taken to 
adopt and test physically motivated interpretations of differ- 
ent functional forms that could be considered when fitting 
these profiles. For example, the lower left panel of Fig- 
ure [T] shows the results of fitting a pure Sersic function to 
the entire surface density profile (including the central star- 
burst component). There is a reasonable fit to the entire pro- 
file with a single Sersic index «, = 3.56, quite different from 
the Sersic index which describes either the pre-starburst or 
starburst light components. Likewise, if we consider the ex- 
cess light to be that light in the real profile above the pre- 
diction of the best-fit Sersic model, we would infer only a 
tiny extra light fraction ~ 1.0%. Although the fit is techni- 
cally worse, with variance = 0.25, the difference is not 
dramatic, and by many observational standards would be con- 
sidered a good fit. Clearly, however, the results do not have 
the same physical meaning in this case - the "extra light" de- 
termined in this manner is no longer a direct tracer of the 

* We fit the observed points in the same manner as the simulations, weight- 
ing each point with the (quadratic) sum of the same intrinsic ~ 0. 1 mag point- 
to-point variance and the observational errors. These errors are however gen- 
erally small (much smaller than the plotted points shown in the observed 
profiles), so weighting the data points equally gives almost identical results. 
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physical starburst component. The differences grow if we 
add a degree of freedom and fit a "core-Sersic" profile, of the 
form / cx r i + {rb/rryil"' exp{-/7„ [(r^ + r^/r,"]'/""} (e.g. 
iGraham et al. 2003), which behaves as a single Sersic pro- 
file outside of r;, and breaks to a power-law of slope 7 within 
rb- The lower right panel of Figure [T] shows the outcome of 
this fit, which is again good in a pure statistical sense, albeit 
worse than our best fit SersicH-extra light fit (cr^ = 0.20 mag). 
However, here the derived parameters are even less intuitively 
related to the known physical decomposition - the best fit Ser- 
sic index is a much steeper = 3.95 and one actually infers 
that the system is a core galaxy, with missing light relative to 
the best-fit Sersic profile. Furthermore, by comparison with 
our results in Paper I, we find that the discrepancy between the 
physical parameters which accurately describe the outer and 
starburst components and those recovered by the pure Sersic 
or core-Sersic profile fits becomes even worse when the mass 
fraction of starburst component is larger. 

The cause of these differences is that the extra light compo- 
nent blends smoothly with the outer pre-starburst light profile. 
By increasing the central surface brightness, the extra light 
component makes the overall profile appear steeper (concave 
up in the /i-r'/"^ projection), owing to the rise at small r. 
However, the cusp itself does not continue to rise steeply in- 
wards (in most cases), so after steepening the best-fit Sersic 
index to fit the outer part of the extra light component, one 
is often forced to infer the existence of a core in the central 
regions. Again, these fits do not directly reflect the physical 
two-component nature of the profiles (rather reflecting some 
combination of the components, with the extra light no longer 
apparent in an excess with respect to the fit but in the higher 
fitted Sersic indices), but they are not terrible matches to the 
light profile. This emphasizes that a physically motivated pro- 
file must be adopted when fitting a parameterized model to the 
data, if one wishes to translate these parameters into robust 
physical properties. Fortunately, there are some indications 
from observations that the Sersic only and core-Sersic fits are 
not physically motivated as a means to decompose the two- 
component nature of the observed systems. 

First, they are technically worse fits, although the differ- 
ence is not large (and in some rare cases, core-Sersic profile 
is a better match to our simulations than the SersicH-extra light 
profile). With photometry accurate to ^ 0.01 mag arcsec"' it is 
possible to robustly distinguish the quality of the fits shown in 
Figure[Tl especially to note the fact that the errors in the Sersic 
and core-Sersic models are a strong function of radius (an in- 
dication of the less ap propriate choice of fittin g function), but 
with the exception of iKormendy et all (12008 ). our data sets 
do not attain such high accuracy. Second, they begin to fail 
at large radii - however, this is where the true nature of the 
Sersic profile of the outer light component is most prominent, 
so any failure at large r should be especially worrisome. Fur- 
thermore, when we examine the kinematics (e.g. ellipticity, 
boxy/diskyness, rotation properties) along the major axis, one 
can often see a transition in these properties where the extra 
light begins to dominate (see § |5]), whereas we would expect 
no such change if the ability to fit a continuous Sersic profile 
were taken to imply that there is only a single physical com- 
ponent constituting the galaxy. Finally, when fitting a core- 
Sersic profile to a system with significant extra light, a large 
missing light fraction (defined as the difference between the 
core-Sersic fit and the inwards extrapolation of the Sersic por- 
tion of the fit) is sometimes seen, even relative to what is typ- 



ically observed in genuine c ore galaxies (e.g. massi ve, boxy, 
slow-rotating ellipticals; see iKormendy et al]|2008h . This is 
true when the dissipational component is large - but Figure[T] 
demonstrates that it is not always the case, so again, care must 
be taken to employ a physically well-motivated decomposi- 
tion and interpretation. 

We do not claim that a pure Sersic or core-Sersic profile is 
never a physically motivated parameterization of the galaxy 
light profile. However, for gas-rich merger remnants, we 
know from our simulations and have good reason to believe 
observationally that there is some excess light component. In 
these cases, which we investigate here, the results of these fits 
are demonstrably less physically intuitive and can be mislead- 
ing. 

We also emphasize that although there are some superficial 
similarities between oi ir adopted p arametric profil e decompo- 
sition and that in e.g. ICdte et alj (|2006) and Fe rrarese et alj 
( I2OO6I) . the two are in detail significantly different and ad- 
dress very different spatial scales and physical properties of 
the galaxies. Typically, the "outer profile" we refer to ex- 
tends to and beyond (in our simulations) the limits of our 
ground-based photometry, corresponding to physical radii of 
^ 20- 100 kpc, and our "inner profile" refers to the resid- 
ual from a central starburst at scales where a significant frac- 
tion of the galaxy mass becomes self-gravitating (see § |6ll, at 
^ 0.5- Ikpc. We stress again that we are not resolving in- 
wards of the central ^ 30-50pc, and our modeling should 
not be extrapolated to within these rad ii without considerable 
care. In contrast, the "outer profile" in lFerrarese et al.l (l2006h 
is based on the HST ACS profiles, which extent to outer radii 
^ 1 kpc, and their "inner profiles" typically dominate the light 
profile at very smafl radii Q.Q\-0.Q2R^ (-- 10-40pc for 
most of their sample). This is more akin to separating our 
"inner" component itself into multiple sub-components - i.e. 
a starburst stellar component that blends (as we have shown) 
relatively smoothly onto the outer, violently relaxed stars and 
an innermost nuclear component. The authors themselves ad- 
dress this, and denote these nuclear excesses as central stellar 
clusters. In Appendix ICl we demonstrate that such clusters 
are clearly distinct, orders-of-magnitude smaller objects than 
the starburst/extra light components we study in this paper 
Such systems may indeed be present (and could be formed 
in the same dissipational starburst which we model): but if so 
they are distinct subsystems sitting on top of the starburst light 
component, which we do not have the ability to model or re- 
solve in our simulations. Therefore, while the two approaches 
may yield complementary constraints, we caution that our re- 
sults are not directly comparable and are specifically designed 
to trace distinct physical structures. 

5. COMPARISON WITH OBSERVATIONS: EXTRA LIGHT IN 
ELLIPTICALS 

We now extend our analysis to the sample of observed el- 
lipticals and merger remnants described in §[3] Figures [3]fT2l 
show surface brig htness profiles o f objec ts in the Virgo el- 
liptical sample of IKormendy et all (l2008h . in order of most 
massive to least massive.^ For each object, we plot the sur- 

' NGC 4382 is typically classified as a core or intermediate galaxy in 
the literature. We show it because it is also sometimes considered a gas- 
rich merger remnant, based on the presence of ripples and shells formed 
by cold tidal material (e.g. Hemauist & Spereel 1992). a high fine structure 
index (Schweizer & Seitzer 1992), young 2Gyr old) stellar populations 
iTrager et al. 2000; McDermid et al. 2006), disky isophotes and high ellip- 
ticity I K ormendv et al..2008.) . and rapid rotation l Emsellem et alj2007i) . Ex- 
cluding it from our sample, however, makes no difference to our conclusions. 
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Fig. 3. — Surface brightness profiles are siiown for cuspy ellipticals in the Virgo cluster. Open circles show the observations, from lKormendv et al.l )2008t) . 
These are the highest-mass cusp or extra light ellipticals in Virgo2J(rvj 2M,). Top: Observed V-band surface brightness profile with our two component best-fit 
model (solid, dashed, and dotted lines show the total, inner/extra light component, and outer/pre-starburst component). The best-fit outer Sersic index, extra 
light fraction, and variance about the fit are shown. Middle: Colored lines show the corresponding surface brightness profiles from the three simulations in our 
library which correspond most closely to the observed system (shown outside to the gravitational softening length, ~ 3()pc). Dashed hne shows the profile of the 
starburst light in the best-matching simulation. The range of oute r Se rsic indices in the simulations (i.e. across sightlines for these objects) and range of starburst 
mass fractions (/si,) which match the observed profile are showitiSl with the variance of the observations about the best-fit simulation (cr^, in magarcsec"-). 
Bottom: Observed disky/boxy-ness (04) and ellipticity profiles, with the median (solid) and 25 — 75% range (shaded) corresponding profile from the best-fitting 
simulations above. Note that these are not fitted for in any sense. Figures |5ll2l show the other cusp ellipticals in the sample, ranked from most to least massive. 

parameters to vary - i.e. we allow limited rescaling in the sur- 
face brightness of the simulated galaxies, but not their radii 
or other properties. Despite the allowed surface brightness 
rescaling, the best-fit simulations almost always have similar 
total luminosities to the observed system, because they must 
have a similar effective radius in order to be a good match. 
Considering ^100 sightlines to each of our simulations (al- 
though, as noted in Paper 1, the observed surface brightness 
profile varies by a small amount sightline-to-sightline), we 
find the best fit to each observed system. 

In the middle panels of Figures l3lfT2l we show the three sim- 
ulations which most closely match the observed light profile. 
For the best-fit simulation, we also show the profile of the 
stars formed in the final, central, merger-driven starburst, as 
described in § H] We show in the figures the outer Sersic in- 
dices fitted to these simulations, along with the typical range 
both across sightlines and across the best-fitting simulations 
(which together give some rough approximation to the range 
of «s which might be observed for these galaxies along differ- 
ent sightlines). We also show the best-fit starburst mass frac- 
tion, along with the range across the best-fitting simulations 
(described below), and the variance of the observed points 
with respect to the best fit'". 

The values shown in Figures Dll2l are based on comparison only to the 
profiles shown, from lKormendv e t al. (2008). In Table[T] the values represent 
the results from all available data sets, including multiple different observa- 



face brightness profile with the best-fit two component model, 
and the corresponding fitted outer Sersic index and extra light 
fraction. It is a reassuring consistency check that our fit- 
ted outer Sersic indices agree well with those estimated in- 
dependently (and with a slightly different methodology) in 
iKormendv et alj (1200 8') - i.e. accounting for the fact that the 
extra light is an independent component (and with sufficient 
dynamic range to well-resolve the outer profile), the results in 
most cases are not highly sensitive to the exact fitting proce- 
dure. We refer to Paper I for the same co mparisons with local 
observed gas-rich merger remnants from lRothberg & JosephI 
(|2004). A complete list of fit parameters and compiled galaxy 
properties is included in Table[T]for the ellipticals and Table|2] 
for the merger remnants. 

We also compare each observed system with our library of 
simulations, in a non-parametric fashion. We do this by al- 
lowing the normalization of the simulated light profile to vary 
(within ±0.5 dex), and quantifying the (variance of the ob- 
served points with respect to the simulated light curve at > 1 
gravitational softening length) of each simulation. We allow 
the normalization to vary because we have a finite number of 
simulations and therefore do not sample a continuum in e.g. 
total brightness, but instead discretely sample at factor ^ 2 in- 
tervals (we do not allow the simulated profiles to vary by more 
than this amount, to avoid an unphysical match to a simula- 
tion with very different total mass). We do not allow any other 
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Fig. 4. — As Figure[3] but in log-log space. The two-component nature of the profiles is somewhat less obvious in this projection, but the dynamic range is 
more clear. We show the systems from the outermost observed radii down to our best simulation resolution limits (~ 30 pc). Over nearly four orders of magnitude 
in radius (and ~ 14 mag in surface brightness), simulations and observed systems agree. We show all of Figures [3ll2l in this projection in Appendix IaI 
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Fig. 6. — The next most massive cusp ellipticals (~ 0.5Mt). 




Fig. 7. — Lower-mass cusp ellipticals (~ 0.2— 0.3M»). Our simulations reproduce the observed outer profiles and kinematic properties of such galaxies, 
do not resolve the stellar cluster nuclei at small radii. The extra light recovered by our two-component fits therefore may be misleading at such low mass. 
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Fig. 10. — Very low-mass cusp ellipticals (~ 0.03 — 0.1 M, ). Our simulations provide less good matches at these luminosities, where dwaif galaxies dominate 
the spheroid population (ellipticals at these masses are very rare). Robustly resolving the extra light in these very small systems probably requires < lOpc spatial 
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Fig. 1 1 . — The lowest-luminosity cusp ellipticals in Virgo (~ 0.01 M,). The comparison with our simulations is similar to Figure [Tol 




In addition, for these simulations we show the isophotal 
shape and ellipticity as a function of major axis radius, com- 
pared to that observed. Note that we do not fit these quantities, 
only the surface brightness profile. We show, for each simu- 
lation, the range across sightlines in these quantities - it is 
clear that these depend much more strongly on sightline than 
the surface brightness profile (this is primarily why we do not 
fit these quantities). In every case, there is a significant frac- 
tion of sightlines with shape and ellipticity profiles roughly 
consistent with those observed, but the simulations highlight 
the range of profile shapes for similar spheroids to those ob- 
served. 

For the intermediate and higher mass Virgo ellipticals, > 
O.IL*, we easily find simulations which provide an excellent 
match to the observed profiles, with variance often less 
than even a multi-component parameterized fit. The fits are 
good over the entire dynamic range from the largest observed 
radii 1 00 kpc) down to our resolution limits 30pc)". At 
radii below our softening limits, the simulation profiles artifi- 
cially flatten; but we show in Appendix |B] that the agreement 
continues down to smaller and smaller radii as we increase 
our resolution. At moderate and large masses, the starburst 
fraction recovered by our two-component fit is usually a good 
match to the physical starburst mass fraction in the best-fitting 
simulations (see §|6]l. 

At the lowest masses L < 0.01 L*, the fitted extra light com- 
ponents tend to be smaller than our simulated starbursts. This 
is at least in part a resolution issue, both in our simulations and 
in the observations (the extra light in the lowest-luminosity 

' ' The dynamic range of the fits is somewhat difficult to discern in Fig- 
ures [3112] owing to the plotting versus r'^"*; we therefore reproduce these 
figures plotting fi versus r in Appendix IaI 



ellipticals is poorly resolved even with HST data). Below 
Mb ^ -18 (i.e. roughly an order of magnitude below ^ L*), 
it is also no longer clear that the fitted extra light components 
are the same physical entities - some of the sharp central fea- 
tures in the profiles may in fact be nuclear stellar clusters, 
modeling of which would require resolving individual star- 
forming complexes in our simulations. (See Appendix|C]for a 
discussion of the differences between extra light c omponents 
studied her ein and the stel lar nuclei studied in e.g. Cot e et af] 
(.2006) and iFerrarese et al]c200 6). Those nuclear components 
are on much smaller scales and bear little resemblance to 
the starburst component we are interested in here.) The two 
lowest-mass "c ompact ellipticals" in our sample, NGC4486b 
and VCC1199 (Bing geh etal][T98l , with R, - 150-200pc, 
are somewhat smaller than any of our simulated merger rem- 
nants. Given our resolution limits and limited sampling of e.g. 
initial disk sizes at these very low masses, it is premature to 
say w hether some different physics (e.g. tidal stripping; Fabe^ 
Il973h is needed, but the systems lie on the fundamental plane 
and their outer profile shapes appear normal. 

Figure [T3] again shows the observed and best-fit simulated 
surface brightness profile s, for a subs e t of con firmed cusp el- 
lipticals in the sample of iLauer et alj (12007 al) . in order of V- 
band magnitude from brightest to faintest. For each simula- 
tion in our library, we have a corresponding to its good- 
ness of fit to the observed profile, and the genuine physical 
starburst mass fraction /sb. We can therefore construct a x^- 
weighted distribution of /sb for each observed system - es- 
sentially, the probability, across a uniform sample of initial 
conditions, that the observed profile was drawn from a simula- 
tion with the given starburst mass fraction. These are shown, 
and compared to the fitted extra light fraction for our two- 
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Fig. 13. — Observed surface brightness profiles of a subset of the confirmed cuspy ellipticals from the sample of lLauer et alj J2007a), with the best-matching 
two component parameterized fit (dashed and dotted lines) and best-fitting simulations (red, orange, and green lines), as in Figure [3] Where multiple sources 
of photometry are available, independent fits to each are shown. The objects are ranked from brightest to faintest in V-band (as shown). Profiles are shown 
over a constant angular scale (top axis; bottom axis shows physical radius in kpc). The corresponding (right) panel for each shows the distribution of physical 
starburst fractions for the simulations which provide a good fit to the observed profile (as described in the text), with the fitted (parameterized) extra light fraction 
(blue dotted line; one for each source of photometry). Where available, red dashed lines show the mass fraction of a secondary (recent) starburst population 
independently estimated from stellar population studies in other works. Our simulation resolution limits do not extend within the central ~ 30 — 50pc, and our 
fits are not intended to describe these radii. 

component models. In general, the fitted extra light fraction 
corresponds well to the characteristic starburst mass fractions 
in simulations which produce similar light profiles. 

There are a small number of observed objects for 
which detailed spectral energy distributions (SEDs) have 
enabled two-component stellar population models to be 
fit, in which there are generally an older, smoother 
distribution and a younger, metal-enriched single burst 
population jTitus et alj 119971; ISchweizer & Seitzeij 119981: 
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'Schweizer 1996; Sc hweizer & Seitzeri 120071; iReichardt et"!!] 
2001; Michard 2006). For these objects, we plot the ob- 
servationally estimated mass fraction (from these studies) in 
the secondary burst population (the vertical red dashed lines), 
which we expect should correspond (roughly) to the starburst 
population in the spheroid-forming merger, if this is a good 
description of the formation history. Although there are only a 
few systems for which this comparison is possible, the agree- 
ment between this estimate and our inferred extra light or star- 
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burst fractions is surprisingly good. Of course, there are a 
number of uncertainties and degeneracies in an attempt to ob- 
servationally decompose stellar populations, but this gives us 
confidence that there is physical meaning to our decomposi- 
tions. 

6. PROPERTIES OF "EXTRA LIGHT" PROFILES 

Having fit both our simulations and observed cuspy ellip- 
ticals to an outer violently relaxed component and an inner 
starburst profile, we now compare these fits as a function of 
galaxy properties. 

6.1. Outer Profiles: Sersic Indices 

Figure [14] plots the outer Sersic indices of our sample of 
simulations as a function of galaxy stellar mass, compared 
with those observed. The lack of any trend is striking - we 
predict that there is no dependence of outer Sersic index on 
galaxy mass for cusp ellipticals. In fact, we have searched our 
entire sample of mergers attempting to find a dependence of 
outer Sersic index on some galaxy property, including merger 
redshift, gas fraction, halo concentration, baryon fraction, and 
the presence or absence of initial bulges, and find no depen- 
dence. There is a weak trend with gas content, but only in the 
sense that systems with extremely high gas content even at 
late merger stages (e.g. > 40% gas at the time of final coales- 
cence) can form or retain massive disks, bringing them closer 
to nj = 1 . There is also a weak dependence on orbital parame- 
ters, but only in the sense of different, extreme orbits changing 
the best fit outer Sersic index by A«j < 1 . The apparent dif- 
ference between our low and high-redshift simulation dis- 
tribut ions is in f act entirely attributable to these effects. Simi- 
larly, |N^S]&^uiillo (2006) find that in simulated collision- 
less (gas-free) disk merger remnants - i.e. systems for which 
the entire profile is by definition part of the "outer," violently 
relaxed component - there is also no significant dependence 
of the Sersic index on mass, effective radius, or merger mass 
ratio. 

The observat i ons appear to confirm this prediction. 
iKormendv etH] ( l2008h see no dependence of outer Sersic in- 
dex on galaxy luminosity (within the extra light/cusp popula- 
tion), and our other data sets support this over a large baseline 
in luminosity and stellar mass (albeit with larger scatter, ow- 
ing primarily to the lower quality of the data). Over more 
than three orders of magnitude in stellar mass, and two orders 
of magnitude in effective radius, the observations and simu- 
lations both show a typical n, ~ 2 - 3 with no dependence of 
outer Sersic index on mass, luminosity, or radius in cusp el- 
lipticals. 

This prediction appears to contradict some previous results 
that argue for a strong dependence of Sersic index on luminos- 
ity ((Graham 2001; Trujillo et al. 2002; Ferrarese et al. 2006) 
or ef fective radius (ICaon et al.l 119931: IPrugniel & SimienI 
Il997h . However, we emphasize that these fits are not di- 
rectly comparable to ours. First, these correlations were 
found considering samples of a broad range of spheroids - 
from dwarf spheroidals through cuspy, rapidly rotating el- 
lipticals through massive, cored, slowly rotating ellipticals. 
Here, we are only arguing that the specific subclass of cuspy, 
true ellipticals formed in gas-rich mergers should have a 
mass-independent Sersic index distribution. Different for- 
mation mechanisms (for example, subsequent dry mergers; 
iHopkins et al.l20 0^h) can systematically change the Sersic in- 
dex, giving rise (via cosmological trends towards more merg- 
ers in higher-mass systems) to mean correlations between Ser- 



sic index and galaxy mass or size. To the extent that the ob- 
served Sersic indices of the cusp population are relatively low 
and do not depend on mass or luminosity, it implies that they 
are generally formed in a small number of major mergers, 
without substantial subsequent re-merging. 

Second, these authors were often fitting the entire galaxy 
light profile to a single Sersic or core-Sersic law, whereas we 
have attempted to decompose the inner and outer galaxy light 
profiles. This two component approach will systematically 
yield different Sersic indices, and in some cases as demon- 
strated in §|4]the difference can be dramatic. 

To check if the effect illustrated in Figure [14] is caused by 
large scatter in our estimates (possibly obscuring an under- 
lying trend), we show in Figure [15] the same Sersic index as 
a function of mass, but using a strict, cleaned sample. We 
include in this subsample only those galaxies for which > 3 
sources of photometry (or photometry in > 3 bands) yield n^ 
values different by < 20%. We exclude the recent merger 
remnants, for which unrelaxed features may introduce addi- 
tional scatter or uncertainties. For the simulations, we plot 
only those with sightline-to-sightline variance less than this 
amount (usually eliminating those with significant tidal or un- 
relaxed features). The distribution and lack of dependence on 
mass is the same as Figure[T4] but the conclusion here is more 
robust. 

Given that, at least for cuspy ellipticals formed in gas-rich 
mergers, there is no strong systematic dependence of n^ on 
other properties or initial conditions, it is acceptable to place 
all such systems on the same footing and consider the overall 
distribution of values. Figure [16] shows this, for the simu- 
lations and our observed samples. Figure [TT] shows the same, 
but restricted to the cleaned subsample of Figure [15] In each 
case, there is rea sonable agreement, within the errors. 

The sample of iRothberg & JosephI (l2004l) shows somewhat 
more n, ^ 1 systems than our simulations or observed cusp 
ellipticals. This effect is only marginally significant, but prob- 
ably arises because a few of the systems in their sample will 
most likely (once they are relaxed) be better classified as SOs 
than ellipticals. There is a tentative suggestion that our pre- 
dictions are shifted to systematically higher Wj than the ob- 
served systems, by An ~ 0.25 or so. At this level, however, 
observational issues in the measurements become important, 
as do the exact orbital parameters used in the simulations 
and the dynamic range oyer wh ich the fit is p erformed (see 
e.g. 'B ovlan-Kolchin et al.ll2005b . For example. iBlanton et alj 
(2005) find that for SDSS light profiles, the observations may 
be biased to underestimate b y ~ 0.2-0.5, an d are sensitive 
to the sky subtraction (see also iLauer et al.ll20b7b) . consistent 
with the offset we see. It is therefore not surprising that the 
agreement is not exact. The important thing is that, with few 
rare exceptions (< 5% of cases), cusp y ellipticals and gas-ric h 
merger remnants have < 4 (as in iKormendv et alj |2008), 
i.e. are concave-down in /i-r'/'* space, without a significant 
dependence of on other galaxy properties. 

6.2. Dissipational ("Extra Light") Mass Fractions 

Figure[T8]compares our estimates of the dissipational mass 
fraction in the observed ellipticals: the directly fitted extra 
light fraction /extra and inferred starburst mass fraction /sb 
from the best-fitting simulations. For clarity we restrict to the 
cleaned sample from Figure [15] Our fitted decomposed extra 
light fraction reliably traces the inferred starburst mass frac- 
tion, with a factor ^ 2 scatter similar to that predicted from 
our simulations (Figure [2]). The starburst fraction /sb itself 
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remnants have characteristic n ~ 2—3, without a strong systematic dependence on mass or other properties. Top: Simulations: color encodes gas fraction, symbol 
encodes orbital pai'ameters (the two "random" orbits are two different common random orbits, one somewhat closer to prograde, the other - "alt" - somewhat 
closer to polar), and filled/unfilled encodes the initial redshift of the simulations (as in plotted key). The plotted ris is the median across ~ 100 sightlines, typical 
sightline-to- sightline difference s are shown as the filled plotted error bar. Bottom: Open colored points are obs erved systems: red stars are the cusps in the Virgo 
elliptical sample of Kormendy et al.| 1200 3) . purple circles the confirmed cusps in the local elliptical sample of Lauer et alj (2007a) (supplemented by the sample 
of Bender et al. ( 1988)), and cyan squares are the gas-rich merger remnants from Rothberg & Joseph 1 2004). Open plotted en'or shows the typical differences in 
derived from different sources of photometry and/or different observed wavelengths. We use this point notation throughout. 




must, in some sense, reflect the cold gas mass available in 
the disks just before the final merger (and we show in Paper 
1 that this is the case) - in this physical sense, our fitted /exua 
and inferred /sb are a robust reflection of the gas content of 
the progenitors. Of course, changing simulation properties 
such as the presence or absence of an initial bulge, the con- 
centration of the progenitor halos and disks, the presence or 
absence of a supermassive black hole, and the treatment of 
star formation and the ISM equation of state can indirectly in- 
fluence /sb by altering how efficiently gas is consumed and/or 
expelled before the final merger, and therefore how much is 
available to participate in the starburst. For a fixed gas mass 
at the time of the final starburst, however, the starburst com- 
ponent mass (and therefore also extra light mass, which traces 
the starburst) is independent of these effects. 

Figure [18] also compares the results of our fitting to in- 
dependent stellar population-based estimates of the starburst 
fraction in observed ellipticals. It is in principle possible, 
by studying the stellar populations in sufficient detail, to 
estimate the mass fraction which formed in a recent, cen- 
tral starburst (as opposed to the more extended quiescent 
star formation history), and this should provide an indepen- 



dent check of our decompositions. Unfortunately, there are 
still a number of degeneracies, and this requires detailed ob- 
servations, but it has been attempted for several of the ob- 
served systems ([ Titus et al. 1997; Schweizer & Seitzer 1998t 
ISchweizer"1996; Sch weizer & Seitzer. .2007; .Reichardt et al] 
2001; Michard 2006). Comparing our estimated /extra or /sb 
with these estimates for the mass fraction in the secondary 
(newly formed/starburst) stellar populations, we find a rea- 
sonable correlation. Although there are only a few objects for 
which sufficiently accurate stellar populations are available to 
allow this comparison, they all suggest that our fitted extra 
light component is indeed a good proxy for the mass fraction 
which was involved in the central, merger-driven starburst. 

Having some confidence that our estimates of /sb are rea- 
sonable. Figure [19] plots the inferred starburst mass fraction 
/sb for the observed systems as a function of stellar mass. 
In the same manner that we have defined a best-fit /sb from 
the best-fit simulations, we can also define a best-fit "initial" 
gas fraction (roughly the gas fraction 1 Gyr before the fi- 
nal merger), and show this as well. We emphasize though 
(for the reasons above regarding the efficiency of pre-merger 
gas consumption and expulsion) that this is a much less ro- 
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Fig. 15. — As Figure [T4l but using a cleaned sample: we include only 
those galaxies for which > 3 sources of photometry (or photometry in > 
3 bands) yield n, values different by < 20%. For the simulations, we plot 
only those with sightline-to-sightline variance less than this amount (usually 
eliminating those with significant tidal or unrelaxed features). Dotted line 
shows the median ris ^ 2 .6. The distribution and lack of dependence on mass 
is the same as Figure [T4l but the results here are more robust. 

bust quantity. In either case there is a clear trend of increasing 
dissipation (increasing fractional mass required in a dissipa- 
tional starburst component) at lower masses. The significance 
of the correlation is unambiguous (> 8 it). We can conve- 
niently approximate the trend in dissipational mass fraction 
as a function of stellar mass with the fitted function 
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Fig. 16. — Distribution of outer Sersic indices in cuspy ellipticals, us- 
ing our two-component decomposition. Solid black line shows the result for 
our entire sample of simulations (each across ~ 100 sightlines). Colored 
lines show the results for the observed samples of Lauer et al. ( 2007a) (top), 
Kormendv et al. (2008) (middle), and Rothberg & Joseph (2004) (bottom), 
with Poisson error bars. iRothberg & Joseph ^2004) include some likely SOs, 
yielding a larger fraction of n = 1 systems. Cuspy ellipticals have a fairly nar- 
row range of n ~ 2.50— 2.75±0.75, in good agreement with gas-rich merger 
simulations. 




with roughly a constant factor ~ 2 intrinsic scatter at each 
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Fig . 17. — As Figure [T6l but for the robust cleaned sample of Figure [Tsl 



Admittedly, the estimation of this trend requires some com- 
parison with our simulations, and one might argue that per- 
haps it is driven by some deficiency in them. However, we can 
repeat this exercise with the empirically fitted extra light com- 
ponent/extra, and show the results in Figure|20l The trend seen 
in /extra(A^*) IS Completely consistent with that in /sb, but with 
a scatter larger by a factor 2 (exactly what we expect, based 
on the predicted and observed scatter in /extra (/sb))- Consider- 
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Fig. 18. — Top: Comparison of our estimated mass in tlie fitted extra ligiit 
component (/extra) versus t he star burst mass fraction in tlie best-fitting sim- 
ulations (/s[,; as in Figures [3|13) . For clarity, we show only observations 
from our "cleaned" sample and where it would be possible for our simula- 
tions to resolve the extra light. In these cases, the two estimates agree well, 
with a factor ~ 2 — 3 scatter in /extra(/sb) (similar to what we expect from 
our simulations; see Figure |2). Middle: /extra versus independent observa- 
tional estimates of the mass fraction formed in a more recent starburst/star 
formation event, from two-component stellar population model fits to the ob- 
served SEDs (Titus et al. 1997; Schweizet 1996; Schweizer & Seitzer 1998, 
|2007|; iReictodt et alji20Qli ; iMichardi2006) Bottom: Same, but comparing 
the stellar population estimates to /ji, from the best-fitting simulations. More 
observations are needed to independently test our estimates, but the stellar 
population data independently suggests our decompositions are reasonable. 

ing just the data in Figure l20l even given its increased scatter, 
the trend of decreasing extra light fraction with mass is signif- 
icant at > 5 (T. We have experimented with alternative, non- 
parametric (albeit less accurate) estimators based on e.g. the 
concentration indices or stellar populations of our simulations 
and observed systems, and obtain a similar answer In short, 
even without reference to our simulations, however we derive 
an estimate of the dissipational component, it is difficult to 
escape the conclusion that it is more prominent in lower-mass 



ellipticals. This confirms a long-standing expectation of the 
merger hypothesis that if spirals are indeed the progenitors 
of ellipticals, more dissipation is required in lower mass sys- 
tems in order to explain their densities and fun damental plane 
correla tions (we examine this in more detail in lHopkins et alJ 
llOOSah . 

This, in fact, should be expected. It is well-established that 
the gas fractions of spirals are strongly decreasing functions 
of mass, at any given redshift. To the extent that these are the 
progenitors of the cusp ellipticals, then, the amount of dis- 
sipation involved in the formation of ellipticals should reflect 
this trend. We therefore compare in Figure[T9]the range of ob- 
served gas fractions of spirals as a function of baryonic mass, 
estimated at z = 0, z = I, and z ^ 2 - 3. The gas fractions fol- 
low, at each redshift, a similar trend to that we find for the 
dissipational fractions of the observed ellipticals. They do, as 
expected for almost any reasonable cosmological history, in- 
crease systematically with redshift, and the typical disk gas 
fractions at z = and z ^ 2 - 3 appear to roughly bracket the 
low and high end of the dispersion in the inferred elliptical 
progenitor gas fractions. 

In other words, the distribution in progenitor gas frac- 
tions implied by the elliptical surface brightness profiles is, 
as a function of mass, exactly what would be predicted if 
one assumes that the progenitors were spirals, and that most 
of the systems were formed by a major merger sometime 
between a redshift of ~ 0-3 Indeed, this is exactly 
what is inferred for the formation times of cusp ellipticals 
from both observations of the early-type or red galaxy mass 
functions (Bundv et al. 2006; Borch et al. 2006: Fontan a et aH 
[2006; Hopkins et al. 2006c, 2008b), from direct stel lar pop- 
ulation synthesis studies (ITrager et al.l 120001: iThomas et alj 
l2005t iGallazzi et al.l l2005h . and by association of elliptical 
galaxy formation with the triggerin g of quasar activity (e.g. 
iHopkins et alll2006aLl2007dL 120081 . 

6.3. Size of the "Extra Light" Component 

The extra light components in simulations and cusp ellipti- 
cals also appear to follow a similar size-mass relation, shown 
in Figure |2T| The correspondence is especially close if we 
consider simulations with initial gas fractions ~ 0.2 - 0.4, 
which tend to be the best analogs to the observed systems. 
This radius is the effective radius of the fitted inner extra light 
component, and is not the radius at which the system appears 
to deviate from the outer Sersic law (but is more physically 
robust). We do not see simulations with extra light effective 
radii ^ lOOpc, corresponding to the smallest extra light com- 
ponents seen in the very low-mass observed systems, but as 
discussed in Appendix |Bl this probably owes to our resolu- 
tion limits. 

In Paper I we show that the size-mass relation is driven 
by the condition that the gas collapsing into the central re- 
gions in the final starburst becomes self-gravitating, i.e. that 
(GMextra/^extra) ~ {GM^,/Re) in tcrms of the extra light mass 
and effective radius. We show that the observations obey a 
similar condition, with (small) scatter and dynamic range sim- 
ilar to that in our simulated mergers. That observed systems 
follow a similar correlation suggests both that we are at least 

Technically the post-merger elliptical mass is not exactly the mass of 
a single initial spiral, but correcting for this amounts to a small horizontal 
shift of the disk and elliptical points relative to one another in the figure (i.e. 
the disk points should be shifted by 0.3 dex to larger masses if all mergers 
are 1 : 1, or 0.1 dex for more likely 3 : 1 mergers), and does not change our 
compaiison (in fact it makes the agreement slightly better). 



0.8 
0.6 
0.4 

0.2 
0.0 



■&» * o 

o o 

o TiWr it 

„ 

o o " -iSr 
^ o 



Cusps in Elliptical Galaxies 
0.8 

0.6 
0.4 



21 



o ooaiD oo 



0.2 
0.0 



.0 Ci 



ID 

T — "T* o o 



'•sWr 



o 



o oocnD OO 



0.1 





Observed Disk 
Gas Fractions: 
z = 

h«H Z = 1 

7 r.: 9-'^ 



9.0 9.5 



10.0 10.5 11.0 
log( M,/M®) 



11.5 12.0 9.0 9.5 



10.0 10.5 11.0 
log( M,/M®) 



11.5 12.0 



0.1 r 



Fig. 19. — Inferred gas content (dissipational/starburst fraction) of elliptical-producing mergers as a function of stellar mass. Initial gas fraction (top) and 
physical final starburst mass fraction (bottom) corresponding to the best-fit simulations to each observed system in the samples of Lauer et al. 1 2007a) (circles) 
and Kormendy et al. (2008) (stars) are shown, with the typical 25-75% allowed range (error bar). Solid line shows the fit to th e data (Equation[3). Colored points 
with error bars indicate the mean (and ±1 cr range in) disk gas fractions at the same stellar mass, at z = ( Bell & de JonH'2001; Ka nnappani2004tiMcGaugli2005t 
blue diamonds, squares, and circles, respectively), z = 1 (Shapley et al. 2005, green squares), and z = 2 (Erb et al. 2006, orange circles). There is a clear trend 
of increasing dissipation required to explain elliptical profiles at lower masses (significant at > Scr), in good agreement with the observed trend in progenitor 
disk gas fractions over the redshift range where cusp ellipticals are formed, and with what is invoked to explain the observed densities and fundamental plane 
correlations of ellipticals (e.g. iKormendv»1989ciHemquist et al.ll993l) . 

by resolution effects over most of the mass range of interest. 

7. IMPACT OF "EXTRA LIGHT" ON GALACTIC STRUCTURE 

7.1. Galaxy Sizes and the Fundamental Plane 

Figure l22lshows how, at fixed mass, the effective radius (of 
the entire elliptical profile) scales with extra light mass. We 
consider three mass bins, below, at, and above ^ . In each, 
we plot Re relative to that expected for the given stellar mass, 
as a function of the fitted extra light fra ction. Spec i fically , 
we determine (7?p(M*)) from the sample of lShen et alj (120031) . 
and take the ratio of the half mass radius of each system (de- 
termined directly from the light profile, or from the fits, it does 
not change the comparison) to that value. Our mass bins are 
small enough, however, that this makes little difference com- 
pared to just e.g. considering in a given bin. There is a 
strong trend: at a given stellar mass, systems with larger extra 
light have systematically smaller R^ (they also have slightly 
larger velocity dispersion a, although the scatter is larger there 
in both simulations and observations). In each case, the sim- 
ulations and observed systems occupy a similar locus. We 
can also construct this plot with the starburst mass fraction /sb 
of the best-fitting simulation as the independent variable, and 
find an even tighter correlation of the same nature. 

This directly implies some structural change (some subtle 
non-homology, albeit not necessarily traditional structural or 
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Fig. 20. — As Figure [19j but showing our empirically fitted /cxna as a 
function of stellar mass. The trend of increasi ng d issipation at lower masses 
is still clear and is consistent with that in Figure [l9l but with an expected extra 
factor ~ 2 — 3 scatter from the scatte r in our purely empirical estimator. Solid 
line shows the best-fit from Figure 1191 dotted lines show the ± 1 cr scatter 
expected based on the scatter in /extra(/sb) (see Figure|2j. Dashed line is a fit 
just to these data (statistically consistent with the solid line, and ruling out no 
dependence of dissipation on mass at > 5 cr confidence). 



roughly capturing the most relevant physics determining the 
scales of extra light, and that we are not being severely biased 
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Fig. 21. — Top: Effective radius of the extra light component {not equiv- 
alent to the radius where it breaks fro m th e outer Sersic fit) as a function 
of extra light mass (points as in Figure [14). Simulations (left) and observa- 
tions {right) are similar, especially if we restrict to simulations with initial 
gas fractions ~ 0.2 - 0.4. Resolution limits (see Appendix |B) prevent us 
from simulating systems with Scxira ^ lOOpc, but this is only important for 
the few very lowest-mass ellipticals (L < 0.01 L,; discussed in §[5). Filled 
diamond is typical sightline-to-sighthne variance in the simulations, open di- 
amond the source-to-source (or band-to-band) scatter in observed profile fits. 
Bottom: Effective velocity dispersion of the extra light component vs. that for 
the whole galaxy. Solid line shows (GMcxtia/^cxtra) = (GMf/Re) - the extra 
light collapses to the point where it is self-gravitating. The observed systems 
{right) follow a trend which agrees well with the simulations (left). 

kinematic non-homology) in the fundamental plane tilt. At 
fixed mass, smaller systems are so because a larger fraction 
of their mass is formed in a central dissipational starburst. 
This dissipational starburst is compact, so even though the 
pre-existing stars are scattered to large radii, the effective ra- 
dius is smaller. 

Given two progenitors of known size and mass, it is 
straightforward to predict the size of the remnant of a dissi- 
pationless merger, simply assuming energy conservation (see 
e.g. [Barnes 1988); in the case of a dissipative merger, we can 
very crudely model the results in Figure |22] by assuming the 
non-extra light stars follow a .Hernquist (1990) profile with 
effective radius Reifsb = 0), and the extr a light is all at r = 0. 
More accurately. ICovington et al.l (l2008h use the impulse ap- 
proximation to estimate the energy loss in the gaseous com- 
ponent, followed by collapse in a self-gravitating starburst. 
This yields a detailed approximation as a function of e.g. ini- 
tial structural and orbital parameters, but if we assume typical 
progenitor disks and parabolic orbits, it reduces to the remark- 
ably simple approximation 



ReiU = 0) 



(4) 



l+(/sb/./0)' 

where /o « 0.25-0.30 and Reifah = 0) is the radius expected 
for a gas-free remnant. We plot this in Figure |22l with the 
scatter seen in the simulations. At all masses, in both simu- 
lated and observed cusp ellipticals, more dissipational ellip- 
ticals are smaller in the manner predicted. In the absence of 
dissipation, the stellar light of observed systems would follow 
a uniform virial relation, but dissipation results in smaller R^ 
at fixed stellar mass, and (given the concentration of mass in 
this central starburst) therefore a higher baryon fraction inside 
Re - i.e. changing the total M/L (total dark matter plus stellar 
ma ss to stellar m ass ratio) within R^ of the stellar light (see, 
e.g. lKormendv& Gebhardt 2001). 

We directly illustrate this homology-breaking in Fig- 
ures I23ll25l by considering the light profiles of systems in 
Figure |22l along the correlation between R^ and /extra at fixed 



stellar mass. We consider the ^ M,^ mass bin (the other two 
give similar results, but have fewer observed systems), and 
identify four regions in the R^ - /extra space along the mean 
correlation, highlighted in Figure |23] We then plot the sim- 
ulated (Figure l24l) and observed (Figure |25] ) light profiles in 
each region. We also plot an = 3 Sersic profile, which pro- 
vides a good fit (in the mean) to most of the outer Sersic pro- 
files in this mass range. 

It is clear in Figure |25] that both simulated and observed 
systems at fixed mass with the largest R,, show profiles close 
to a pure Sersic law, with little evidence for any extra light 
component in their centers (indeed, they have f^\, < 0.03). 
The observed systems in this regime are still cusps, but the 
cusps tend to be prominent at very small radii and (in several 
cases) somewhat shallow, and contribute negligibly to the stel- 
lar mass. However, moving to smaller R,,, deviations from a 
Sersic law at r ^R^ become more prominent. That is not to 
say that these deviations are universal (that the extra light al- 
ways takes the same shape/form), but there are increasingly 
prominent central light concentrations. If the systems were 
strictly homologous, there should be no differences in Fig- 
ure |25] Since we have scaled each system at its Rg, they 
should be identical - instead, it appears that a central light 
concentration, formed in our simulations via gas dissipation, 
drives an important non-homology. 

We discuss this in Hopkins et al. ( l2008ah . where we show 
that this is sufficient to explain the tilt of the (stellar mass) 
fundamental plane. In that work, we demonstrate that the fit- 
ting of the central component in the light profile can be used 
as a direct observational test for the role of dissipation in the 
fundamental plane. Given a larger central light component, 
the effective radius of the stellar light profile must be smaller. 
Because the central regions of the galaxy are the most baryon- 
dominated, moving the effective radius inwards results in a 
different (larger) ratio of stellar to dynamical mass, i.e. driv- 
ing tilt in the fundamental plane. We note that although this is 
technically non-homology (i.e. the profiles are not perfectly 
self-similar), it does not drive tilt in the sense of traditional 
structural or kinematic non-homology; i.e. the practical "ho- 
mology assumption" that the true mass enclosed in R,, is pro- 
portional to the dynamical mass estimator a^R^/G remains 
true. Rather, the non-homology induced is a subtle effect that 
provides a tracer of dissipation, which changes the physical 
ratio of baryonic to dark matter within R^. 

7.2. Galaxy Shapes and Kinematics 

Figures |26] shows how the shapes and global kinematic 
properties vary with extra light. We plot the rotation, isopho- 
tal shapes, and ellipticity of our simulations and the observed 
systems as a function of the fitted extra light fraction. Here, 
we consider the global properties of the galaxy - namely 
the rotation (V/cr)* measured within /?e, and boxy/disky- 
ness 100fl4/fl, and ellipticity e measured for the half-mass 
projected (Rg) isophotal contour, for each of ^ 100 lines- 
of-sight to the remnant uniformly sampling the unit sphere 
(i.e. representing the distribution across random viewing an- 
gles). The details of the fitting for our simulations are de- 
scribed in ICox et al.l (2006b), and f o r the observa t ions i n 
iBender et al.' H9S1); Bender (1988'); 'Bender et alj (Il992h ; 
Faberet al. (1989); Simien & Prugniel (1997, 200% We de- 
fine rotation in the standard manner, relative to that of an 
oblate isotropic rotator, with the parameter (V/a)*, defined 
using the maximum circular velocity Vc velocity dispersion 
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Fig. 22. — Effective radius Re relative to the median value for all ellipticals of the same stellar mass, as a function of our fitted extra light fractions (the 
empirical tracer of the dissipational/starburst mass fraction). We compare simulated gas-rich merger remnants {top) with observed cusp ellipticals and gas-rich 
merger remnants (bottom), with points as in Figure[T4] We show this in three bins of stellar mass (relative to M, fa lO" M(y,, or M*, = ~21) . Solid (dotted) lines 
show the mean (±1 cr) correlation, following the analytic solution for dissipational mergers and fits to our simulation in lCovington et al] (2008). Simulations 
and observations exhibit the same behavior: systems with smaller Re at fixed mass have systematically larger extra light fractions {> 6 a significance in the 
observations). This implies that, at fixed mass, systems are driven along the fundamental plane by the relative amount of dissipation involved in their formation. 
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Fig. 23. — Center panels of Figure l22l with four regions along the cor- 
relation between effective radius and extra light fraction at fixed mass high- 
lighted The surface brightness profiles in these regions are shown in Fig- 
uresl24l&|25] 

within Re, and ellipticity as (lKormendvll982l) 
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(V/ar = (V/a)/^e/(l-e). 



(5) 



We exclude the coplanar merger simulations from our com- 
parisons here: those simulations are idealized perfectly copla- 
nar prograde orbits, and as such produce pathological orbit 
and phase structure (we do, however, include some more rep- 
resentative orbits in that are not far from coplanar). 

In each sense, the simulations and observations occupy a 
similar locus. While there are no real correlations, there some 
broad differences between the distributions for systems with 



Fig. 24. — Light profiles of simulated systems with different Rc at fixed 
stellar mass. We plot all pr ofil es from our simulations (colors denote simula- 
tion gas fractions as Figure fTil from each of the labe led regions in the R^ -fsb 
space for the ~ M* galaxy mass bin in Figure |23| in order from largest Rc 
(smallest fsb) to smallest R^ (largest fsb) (left to right, top to bottom). Similar 
results are obtained for the other mass ranges. The dashed black line shows 
a constant Sersic profile, the same in each panel, for comparison. There is 
a substantial systematic difference: smaller Rc systems at fixed stellar mass 
have more prominent central mass concentrations, driven by dissipation in 
our simulations. 

significant or insignificant extra light components. Systems 
with small extra light fractions (/sb < 0.1) have a wide range 
of (y/cr)*, 04/0, and e - they range from relatively slow ro- 
tation (a lbeit not as slow as m any core elliptical slow rota- 
tors; see lEmsellem et af]|2007l) to rotationally supported ob- 
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Fig. 25. — As Figure l25l but for the obse rved surface brightness profiles. 
Color denotes the observed sample as Figure [141 



jects (y/cr)* > 1, with both boxy and disky isophotal shapes 
(-2 < 100 04/(3 < 2), and a range in ellipticity from fairly 
spherical to highly flattened e ~ 0.5. Despite in some in- 
stances having e.g. slow rotation properties, these systems 
do not follow other trends of massive, cored slow rotators 
- they usually are not simultaneously slowly rotating, boxy, 
and round. This is similar to several observed systems, which 
generally show evidence for formation in a gas-rich merger 
but owing to sightline effects may be deviant in one of these 
properties. 

At high extra light fractions, however, both our simulated 
systems and those observed tend to be more rapid rotators 
(fractionally more populating {V /a)* > 0.4), with less boxy 
isophotal shapes (typical -0.5 < 100a4/a < 2), and slightly 
rounder ellipticities (e < 0.4; note they are still more ellipti- 
cal than core galaxies, in agreement with observed trends). 
A more detailed analysis of the role of central mass con- 
centrations and dissipation in shaping the orbi tal structure of 
as-ric h merger remnants will be the focus of iHoffman et alj 
2007h . Briefly, however, these differences are physically 
expe cted, and follow what has been seen in earlier work 
(e.g. Barnes & Hernquist 1996; Naab et al. '2006a; Cox et al.' 
2006b; Onorbeetal.'2006; Jesseit et al. 2007; Burkertetal. 
20071) : the central dissipational mass concentration in these 



systems is highly concentrated, and acts effectively like a 
point mass at the center of the galaxy to much of the mate- 
rial at Rg. When this becomes a large fraction of the galaxy 
mass (/sb > 0.1), the potential becomes more spherical, 
which eliminates some of the phase space d ensity that might 
otherwise support boxy and radial orbits ( | Binnev & Petrou 



1985| ;lHernauist et alJ1993l;lNaab et al.l2006al;lHoffman et al 



20071) . and the central cusp can itself disrupt box orbits as they 
pass near the center (Gerhard & Binney 1985). This results in 
a larger fraction of rotationally supported tube orbits, diskier 
isophotal shapes, and (directly) a rounder overall system. 

If we consider simulations with sufficiently large dissipa- 
tional fractions (/sb 0.1, similar to what we see in the ob- 
served systems), then, our gas-rich merger remnant simu- 
lations match well each of the distributions of e, 04/0 and 
(V/cr)* for cusp y elliptica l s. This is demonstrated in much 
greater detail in ICox et al.l (l2006bh . Here, we highlight that 
there is a significant difference in many of these properties 
between cusp and core ellipt icals, as has bee n well estab- 
lished in previous studies (e.g. lBender et aljT 989; Fab er et al] 
Il997h . It is clear that our gas-rich merger simulations provide 
a good match to the observed distributions of (V/cr)*, 04/0, 



and e in cuspy ellipticals, as we would expect, but they are 
not reproducing the more slowly rotating, boxy, round distri- 
butions characteristic of core ellipticals. This should also be 
borne in mind considering Figure |26l the lack of real corre- 
lations with extra light fraction appears to be true within the 
cuspy population; if core ellipticals are altered through subse- 
quent spher oid-spheroid "dry" mergers, fliese properties can 
be mo dified (iNaab et al.l2006bl:ICox et al.l2008tlBurkert et al.1 
l2007h . and it is expected that these processes will bring the 
predicted distributions into better agreement with those ob- 
served in the cored population. 

7.3. Nuclear Black Hole Masses 

Figure |27] shows how the BH masses scale with 
extra light mass, both globally and at fixed stellar 
mass. We plot all of our simulations which in- 
clude central BHs, and include the cusp ellipticals in 
our sample for which direct kinematic or maser mea- 
surements of a central BH h a ve been possible (com- 
piled from [M agorrian et al. 1998; Merritt & Ferrarese 200 ij; 
Tremain e et al. 2002; Marconi & Hunt 2003 ; Hiiring & Rixl 
20041; lA ller & Richstone 2007). For more details, we refer 
to iHopkins et al.i (i2007b) . There is a correlation in simula- 
tions and observations between BH mass and extra light mass. 
However, this appears to be largely driven by the correlation 
between BH mass and total spheroid stellar mass, which has 
smaller dispersion and a weaker residual dependence on e.g. 
gas fraction or orbital properties. This is expected: if BH 
mass is actually sensitive to the depth of the central poten- 
tial and spheroid binding energy, as argued from the nature 
of the fundamental plane-like c orrelation for BH mass es and 
host properties demonstra ted in iHopkins et al.l (l2007allbl) and 
lAller & Richston j (l2007h . then this is better correlated with 
the total stellar mass setting the potential than the few percent 
of the mass in the extra light component. 

We might expect though, that at fixed total stellar mass, sys- 
tems with larger extra light components, since this formed 
a more compact remnant, would have somewhat deeper po- 
tentials. Given the black hole fundamental plane (BHFP) 
in terms of stellar mass and (i.e. scaling of BH mass 
with R e or, equivalent ly, a at fixed galaxy stellar mass, as 
in Ho pkins et al.l (l2007albi) '). and the scaling of R^ with extra 
light fraction at fixed mass seen in Figure|22] we can estimate 
the dependence of BH mass on extra light fraction at fixed 
stellar mass. Figure [27] shows the residual of BH mass (i.e. 
BH mass relative to that expected at each stellar mass) as a 
function of extra light fraction, compared to this expectation. 
There is considerable scatter, but the simulations and few ob- 
served systems broadly follow the expected trend - namely 
that the cusp mass does not primarily drive Mbh, but that at 
fixed stellar mass, increasing dissipation builds more compact 
remnants which have higher binding energies and therefore, 
via the BHFP, larger BH masses. 

8. IMPACT OF "EXTRA LIGHT" ON STELLAR POPULATIONS 

If a significant fraction of gas dissipates to produce a cen- 
tral mass concentration, and this forms stars in a rapid central 
starburst, we should expect this to leave an imprint on stellar 
popula tion gradients in the remnant (e.g. lMihos & Hernquist! 
Il994bl). 

Figure |28] presents an illustrative example of these ef- 
fects. We consider the elliptical NGC 3377, for which de- 
tailed stellar population gradi ents have been measured by 
ISanchez-Blazquez et al.l (l2007b . The system is a clear cusp 
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Fig. 26. — Top: Correlation between fitted extra light fraction (in both simulations and observed cuspy ellipticals, plotted as in Figure [14) and global 
kinematic properties of the galaxy: rotation (V/ct)*. boxy/disky-ness 100a4/a, and ellipticity e. Median values across sightlines are plotted. Black points with 
eiTor bars show the typical ~ Icr sightline-to-sightline dispersion in each quantity from our simulations. Bottom: Histograms show the distribution in each 
property, for observed systems (thick black line), dissipational simulations which have similar moderate extra light fractions (blue shaded), and gas-poor (nearly 
dissipationless) simulations (/sb < 0.1; red dotted). The simulation distributions uniformly sample each simulation in solid angle over ~ 100 lines-of-sight 
(including the large sighthne-to-sightline dispersion in the top panels). Broadly, cuspy ellipticals and gas-rich merger remnants occupy a wider range in these 
parameters at low /sb, but are uniformly rapid rotators, with slightly diskier isophotal shapes, and slightly rounder ellipticities at large /jb ^0.1. For more details, 
see Cox et al. l 2006b). 



galaxy and rapi d rotator with disky isop hotal shapes, and is 
included in our iKormend v et al.' ('2008') sample (Figure |5]l. 
From our comparison in Figure |5j we select the ^10 simula- 
tions with the most similar surface brightness profiles. These 
generally fall within the range of those plotted in the figure, 
and all provide reasonably good fits to the observed surface 
brightness profile. For each simulation, we then extract the 
stellar populations as a function of radius, viewed from the 
same projection as that which provides the most similar sur- 
face brightness profile to NGC 3377 (although the sightline- 
to-sightline variation is weak). 

We model the emission from each star particle treating it as 
a single stellar population with the formation time and metal- 
licity determined self-consistently from the star-forming gas 
in our sim ulations, and convolving wi th the stellar p opulation 
model s of iBruzual & Chariot (l2003h . assuming a IChabriej 
( 12003b initial mass function. The observed stellar population 
parameters are effectively light-weighted in the optical SED, 
so as a rough proxy for this we extract, in each major-axis ra- 
dial annulus, the B-band light weighted stellar population age 



and metallicity from the simulation. The g-z color is deter- 
mined directly in each annulus. 

Of course, doing this requires that we pick a definite time 
after the merger to view a simulation. We therefore show re- 
sults for three representative times. First, just ^ 10*^ yr after 
the final coalescence of the two galactic nuclei, when the ob- 
ject would likely be classified as a recent merger remnant. 
Second, 3 Gyr after the final merger, which is representa- 
tive of younger ^ cusp ellipticals. Specifically, we choose 
this time because the total light-weighted age, integrated over 
the whole galaxy, at this point matches that observationally 
inferred for NGC 3377 - i.e. when our comparisons are most 
appropriate. Third, we show the remnant after ^ lOGyr, i.e. 
having evolved in isolation for nearly a Hubble time, com- 
parable to the oldest observed ellipticals. At f = 3 Gyr, we 
directl y plot the observed age , metallicity, and color profiles 
from Sanchez-Blazquez et al.l ( 12007 ). At the other times, we 
make the lowest-order reasonable corrections to highlight rel- 
ative evolution in the shapes of the profiles: we add or sub- 
tract the appropriate age difference uniformly from the age 
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Fig. 27. — Top: BH mass as a function of fitted extra ligtit mass, in our sim- 
ulations (left) and observed systems that have direct BH mass measurements 
(right; only a small fraction of our sample have such BH masses measured). 
Bottom: Residual congelation of BH mass with extra light mass at fixed stellar 
mass (i.e. M^w/ {M-g]\(Mt)) versus Mcxua/M,). Solid line is the expectation 
if systems obey the correlation between /si, and /?c from Figure |22] and the 
black hole fundamental plane (Hopkins et al. 2007b). There is a significant 
Mbh ~ Mcxtra Correlation, but it largely reflects the Mbh ~ correlation or 
BHFP. At fixed stellar mass, the (weak) residual trend comes from the BHFP 
as increased dissipation builds a more compact remnant with deeper poten- 
tial, and therefore larger BH. 



profile, and likewise add or subtract the mean color difference 
expected for a single stellar population of the mean observed 
age and luminosity. Again, these leave the profile shapes un- 
changed - we merely shift them by the expected mean to com- 
pare with the shape changes predicted by our simulations. 

At early times (f = 0.1 Gyr), the gradients are strong. The 
stars formed in the central starburst are very young, and 
thus dominate the central light. Unsurprisingly, then, the 
age at the center of the galaxy is approximately just the 
time since the merger, ~ 0.1 Gyr This rises to ^ 3-5Gyr 
in the outermost regions, representative of the ages of the 
stars that were forming before the merger proper began. 
The central population is also the most metal rich, pro- 
ducing a similar strong metallicity gradient. The young, 
blue stellar populations of the center result in a strong color 
gradient with a blue core, common among yo ung merger 
remnants (e.g. [S chweizer & Seitzer 1992; Schw eizerlll996t 
iRothbe rg & Jos epb2004j: jYamauchi & Gotoii2005]) . 

The gradients in these quantities are most pronounced be- 
tween ^ Q.Q5-Q.lRe and ^ R/. at much smaller or larger 
radii they tend to flatten. This intermediate region is the tran- 
sition between dominance of the pre-merger stellar popula- 
tions and central cusp - the gradients are primarily driven by 
the mean difference in stellar population parameters between 
the central starburst and outer old stellar components. This, in 
principle, allows arbitrarily strong gradients at the radii typ- 
ically observed, but at larger radii, for example, violent re- 
laxation mixes the old stars, washing out initial gradients and 
leading to a flattening in the total gradient. As predicted, these 
trends are typical in observed stellar populatio n gradie nts (e.g. 
Mehlert et al.l2003l:ISanchez-Blazquez et al.l2007:,Reda et all 



2007h . 



By ~ 3 Gyr after the merger, the age and metallicity gradi- 
ents have weakened slightly. The apparent age gradient, quan- 
tified as e.g. dlog(age)/dlog(r) is much weaker at this time, 
but largely for artificial reasons - adding a uniform ^ 3 Gyr to 
the age of the system accounts for most of this, since it makes 



the difference in log (age) smaller. In terms of d(age) / d log (r), 
the effects are more subtle, comparable to what is seen in the 
metallicity gradient. The gradients do, at this level, weaken 
slightly. This is because the now older central stellar popula- 
tions have a lower L/M, more comparable to the old stars with 
which they are mixed. There is therefore slightly more mixing 
between the pre-merger and starburst populations. This can be 
seen in the metallicity gradients, which we would expect (to 
lowest order) to remain constant with time. Still, the effect is 
clearly second-order These trends continue with time, as can 
be seen at the time ^ 10 Gyr after the merger. 

The color gradients evolve significantly with time, however 
The inner regions are younger, which at early times (when 
they are very young, < 0.5 Gyr) typically results in blue cores 
and a color gradient which becomes redder at larger radii 
(note, however, that at times very close to the merger, dust 
can reverse this trend). However, they are also more metal- 
rich, which tends to make them redder After ^ 1 Gyr, the 
effects of the age difference are much less prominent, and 
the color gradient becomes dominated by the metallicity gra- 
dient, resulting in red cores, and reversing the sense of the 
gradient. At even later times, the strength of this gradient 
(towards redder central regions) becomes stronger, although 
it remains relatively weak. We therefore expect that most 
ellipticals, even those which are relatively young, wil l have 
weak color gradients with red cores , as observ ed (Faberet alJ 
1989; Bender et al. 1993 HTrager eTa l. 2000; Cote et al.ll2006l: 
IFerrarese et al.. ,2006) . This should caution strongly against 
inferring too much from an observed color gradient, since the 
physical meanings and sense of the color gradients in typi- 
cal simulations not only quantitatively depend on time and 
relative age, metallicity, and abundance gradients, but in fact 
usually reverse their meaning and behavior with time. 

Because the star formation history depends on spatial loca- 
tion, mergers and dissipation can also induce gradients in e.g. 
the chemical abundance patterns and a-enrichment of ellipti- 
cals. However, we do not show these explicitly for two rea- 
sons. First, our simulation code does not currently track sepa- 
rate enrichment by different species, so our estimates of such 
are based on crude analytic estimates (taking e.g. the star for- 
mation history in radial bins and estimating the a-enrichment 
based on a closed-box model). Second, and more important, 
the effects are much more sensitive to the initial conditions 
and cosmological merger history. For example, one can imag- 
ine a situation where stars form over a fairly extended time pe- 
riod in the initial disks, leading to relatively high [Fe/Mg] at 
large radii. The merger then forms a rapid starburst in the cen- 
ter, which will preferentially be a-enriched but not have time 
to self-enrich in heavier elements, leading to a decreasing a- 
enrichment gradient. However, to the extent that most of the 
metal content of the central starburst often comes from the 
previous history of enrichment of that gas in the pre-merger 
disks (i.e. gas that enriched over long timescales), the addi- 
tional effect of the starburst will be quite weak, leading to 
gradients of only ~ 0.1 in [Fe/Mg] over ~ 0.01 - lORe. 

Furthermore, one can also imagine a scenario in which the 
progenitor disks rapidly exhaust their gas in a short timescale 
leading to high a-enrichment, retaining a relatively poor gas 
reservoir for the final merger. Because this gas is retained 
over long time periods, it enriches in Fe-elements from the 
older disk stellar populations, and after forming stars in the 
central starburst creates a gradient with the opposite sense (i.e. 
decreasing a-enrichment with radius). Again, our estimates 
suggest the effect will be weak. 



Cusps in Elliptical Galaxies 



27 




0.01 



0.10 



1.00 



0.01 



0.10 



1.00 



0.01 



0.10 



1.00 



r/R„ 



r/R„ 



r/R„ 



Fig. 28. — Comparison of the observed stellar population and color gradients (points) of NGC 3377 with those in the ~ 10 simulations (lines) which most 
closely match its surface brightness profile (see Figure |5). The B-band light weighted parameters are plotted at three different times after the merger (labeled), 
with At = 3Gyr corresponding roughly to the observed mean stellar population ages in NGC 3377 (points there are plotted exactly as observed). At A? = 0.1 and 
lOGyr, in contrast, we have added the mean age offsets (— 2.9Gyr and +7 Gyr), and shifted the colors by the mean difference expected for passive evolution of a 
single stellar population with solar metalhcity and the quoted age; this highlights the change in shape of the simulation gradients, as well as their normalizations. 
The simulated age and metallicity gradients fade slightly with time, as the difference in stellar populations becomes less prominent. At all times, however, 
the gradients are primarily driven by the mean difference in age and metallicity between the compact central starburst populations and the violently relaxed 
pre-merger populations. At early times, young, starburst populations often yield blue cores and a rising color gradient; at late times, the age effect weakens, and 
the enhanced metallicity of the starburst population dominates, reversing the color gradient (yielding red cores). 



In either case, because the magnitude of the induced a- 
enrichment gradien ts are comparable to those observed in 
spiral galaxies (e.g. Moorthv & Holtzmanll2006t iGanda et alj 
120071) . the final prediction is much more sensitive to the ini- 
tial conditions (unlike e.g. the induced total metallicity gra- 
dients, which are typically stronger than those in comparable 
disks). The mean or total a-enrichment, on the other hand, 
depends of course primarily on the timescale for the forma- 
tion of most of the stellar mass - which in almost all the cases 
of interest here means the pre-merger disks, which we are not 
attempting to model. Further complicating matters, cosmo- 
logical infall or accretion of new gas and minor mergers will 
play a (fractionally) larger role. For these reasons, a more 
comprehensive set of models is needed to say anything quan- 
titative about a-enrichment gradients, but it is nevertheless 
reassuring that the magnitude and sense (both positive and 
negative) of these possible induced gradient s are in fact qual- 
itatively similar to what is o bserved (e.g. iReda et alJ 120071; 
[Sanchez-Blazq uez et al]|2007h . 

Figure|28]is instructive, but it is difficult to compare with all 
systems. We therefore globally compare gradient strengths 
by fitting the mean gradient in metallicity and age over the 
range ^ Q.\-\Re in each of our simulatio ns, comparable to 
the range typically used observationally (e.g. lReda et al.l20O7l 
and references therein). We consider the results extracted at a 
series of times after the merger similar to those in Figure [28] 

Figure |29] shows how these gradient strengths compare 



to those in observed ellipticals (from iMehlert et alj 120031 : 
iReda et al.l2007t ISanchez-Blazquez et alj2007h . We consider 
gradient strength as a function of extra light fraction, velocity 
dispersion, central stellar population age (averaged within an 
aperture of /?e/8, again comparable to typical observations), 
and central metallicity. There is a noticeable correlation be- 
tween gradient strength and the excess light fraction; we ex- 
pect this, since as noted above this central light drives most 
of the gradients. This also gives rise to a correlation between 
e.g. metallicity gradient and central metallicity, age gradient 
and central age, and different gradients themselves (stronger 
metallicity gradients tend to accompany stronger age gradi- 
ents). There is not a strong dependence of gradient on mass, 
however (what there is is mostly driven by the mean depen- 
dence of extra light fraction on stellar mass). In each case, the 
observations occupy a similar locus to the simulations. 

The re have been some claims that the metallici- 
ties (|Naab & Ostriked l2007h and metallicity gradients 
(iForbes et alJl2005l) in ellipticals are too large/strong for them 
to be formed from mergers of present-day spirals. However, 
these arguments fail to account for the role of dissipation in 
increasing the central metallicity and giving rise to strong 
metallicity gradients in merger remnants. The left panels of 
Figure |29] demonstrate that, for moderate gas fractions, our 
simulations occupy a similar locus in both central metallicity 
and metallicity gradients to observed ellipticals. The metal- 
licities typically quoted in stellar population studies are cen- 
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Fig. 29. — Top Left: CoiTelation between metallicity gradient (d[Z/H]/dlogr) and extra light fraction (/extra) in our simulations (points as in Figure[T4j. 
Remaining Panels: Various correlations between metallicity gradients (d[Z/H]/dlogr), age gradients (dlog(Age)/dlogr), central velocity dispersion a, central 
light-weighted stellar population age (Age^,), and central light-weighted stellar population metallicity ([Z/H]^,). For clarity, the simulations are plotted in these 
panels as light grey points, and compared with observed ellipticals (magenta points; triangles, squares, and diamonds are from the samples of Mehlert e t al.l 
jT003); Reda et al. (2007); Sanchez-Blazquez et al. (2007), respectively). The correlation between gradient strength and the excess Ught fraction (top left) reflects 
the fact that the central light drives most of the gradients. This gives rise to the associated correlations shown in other panels. In each case, the observations 
occupy a similar locus to the simulations. 



tral metallicities, measured within ~ Re/^ (or within a cen- 
tral fiber in automated surveys such as SDSS). It therefore re- 
quires only a small amount of material to rapidly self-enrich 
in the central regions to explain these metallicities - our simu- 
lations even with ~ 5% excess light fractions are able to do so 
(reaching [Z/}\\ - 0.5 and d[Z/H]/dlog(r) - -0.6). Since 
strong central metallicities often accompany strong metallic- 
ity gradients, the actual total metal mass of the galaxy (esti- 
mated from the combination of the surface density profile and 
metallicity profile) can easily be a factor ^ 3 smaller than the 
metal mass that would be inferred if the system had a constant 
metallicity equal to the quoted central values. 

We emphasize that the pre-merger stars in our simulations 
are initialized to have a low uniform metallicity with no ini- 
tial gradients (so that we can be sure any gradients we see are 
the result of the simulation dynamics, not artifacts of our ini- 
tial conditions). Since more realistic initial disks will already 
have high metallicities ^ (similar, in fac t, to ellipticals of 
the same mass; see e.g. iGallazzi et al.ll2005h . and metallicity 
gradients in the same sense as those produced in our simula- 
tions, the metallicities and metallicity gradients in our simu- 
lations are only lower bounds to the true values that can arise 
in dissipational mergers. Furthermore, recycling of metal-rich 
gas by stellar evolution (wh ich we do not model here, but see 
e.g. lCiotti & Ostrikei1l2007h will subsequently enrich the sys- 
tem. 

Since the behavior of the color gradients is somewhat more 
complex, we try and reduce it to the key qualitative element. 
Figure l30l shows the distribution in color difference at 
and at Re, a rough proxy for the color gradients, at various 
times after the merger Observed distributions from samples 
of cusp ellipticals, with typical ages 3 - 8 Gyr, are also 
shown. The transition from initial blue cores shortly after the 



merger to red cores in relaxed ellipticals is obvious, as is the 
weakness of the color differences. The simulations and ob- 
servations trace broadly the same distribution for post-merger 
times of ~ 3 - 10 Gyr We find a similar agreement compar- 
ing to the V -I color gradients in Lauer et al. (2005), and the 
V -R and V - / color gradients in B ender et al.l (1198 8 ) . 

The weak color gradients in observed ellipticals are often 
used to justify the assumption of a constant stellar mass-to- 
light ratio as a function of radius, an assumption we have used 
in comparing our simulations and observed systems. In Fig- 
ure[3Tl we test this assumption directly in our simulations. We 
consider two ^ galaxy mergers, one relatively gas poor (a 
case which happens to provide a reasonable match to many of 
the observed ellipticals), and another otherwise equal but gas 
rich merger, both of which have typical color and stellar popu- 
lation gradients. We directly calculate the stellar mass-to-light 
ratios in narrow major-axis radial annuli in both /T-band and 
Z?-band, at different times after the peak of the merger-induced 
starburst. 

The M /L gradients reflect what is expected from the color 
gradients, and supports our assumption of nearly constant 
M/L with radius. At early times ^ IGyr, the central stel- 
lar populations are bright (associated with blue cores), and in 
B-band, nuclear dust obscuration makes the central M/L sen- 
sitive to viewing angle. This is a relevant for the recent merger 
remnants studied in Paper I (although they are observed in 
the K-bwd), and we discuss the implications for the youngest 
systems (with very large apparent /extra ^ 0.5) therein. By 
~ 2 Gyr after the merger, however, the remnants have typical 
weak, red cores and M/L in optical bands such as B and V 
is usually constant as a function of radius to within < 20%. 
At these ages, the parameters derived by fitting to optical pro- 
files are identical to those obtained directly fitting the stellar 
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Fig. 30. — Color gradients (here the difference in mean color within Rf/S 
and Rf) in simulations (black lines) at different times after the merger (la- 
beled). We compare to the distribution for observed cusp elliptical samples 
(for the Cote et al. 1 2006) objects, we show the distribution for ellipticals 
< 6 Gyr [middle^ and > 6 Gyr [bottom] old). Positive values indicate blue 
cores. For ages typical of cusp ellipticals, the remnants have relaxed to where 
metalHcity gradients dominate the overall color gradient and yield red cores, 
albeit with generally weak gradients, with distributions similar to observed 
ellipticals of comparable ages. 
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Fig. 31. — Mass-to-light ratio {M/L) as a function of radius and time 
elapsed since the end of the merger-induced starburst, in S'-band (top) and 
B-band {bottom), for relatively low gas fraction (left) and high gas fraction 
{right), otherwise equal "typical" (~ L,) merger simulations. Shaded range 
shows the ~ Ict range in M/L across different sightlines. At early times 
(< 1 Gyr), young central stellar populati ons lead to a significant gradient (re- 
lated to the strong "blue cores" in Figure [30l . By > 3 Gyr (relevant for almost 
all the ellipticals in our sample), there is no significant dependence in optical 
bands, and only weak dependence in the near infrared (owing to the remain- 
ing metallicity gradients; for detailed discussion see Paper I). 
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Fig. 32. — Correlation between metallicity gradient strength and effective 
radius (relative to the mean at each stellar mass), at fixed stellar mass (in the 
mass bins labeled). Our simulations ai'e show n (p oint style as in Figure [14) 
along with the observed systems from Figure [29] (black circles and squares 
from Sanchez-Blazquez et al. (2007) and Red a et al.. (20 07). respectively). In 
the simulations, increased dissipation yields more compact remnants, with 
stronger stellar population gradients (unfortunately the color and age gradi- 
ents are strongly sensitive to subsequent time evolution effects, so are less ro- 
bust tracers of dissipation), predicting a correlation of the nature shown. Ob- 
servations tentatively suggest the same, but more are needed to confirm this 
prediction of different degrees of dissipation along the fundamental plane. 



mass profiles. The gradients in age (younger at the center) and 
metallicity (more metal rich at the center) have opposite ef- 
fects in the optical, yielding negligible dependence of M/L on 
radius (in K band, the effects do not cancel; however among 
our observed samples only the recent merger remnants from 
iRothberg & JosephI (l2004h are observed in /T-band, so we re- 
fer to Paper I for more details). Reassuringly, this is similar 
to a purely empirica l estim ate based on the color-dependent 
M/L from .Bell et alj (l2003i) and the observed color gradients. 
All of the ellipticals we study have estimated stellar popula- 
tion ages older than this threshold (and weak color gradients) 
and are observed in these optical bands, so our results should 
be robust to stellar population effects. 

Given the dependence of both gradient strength and effec- 
tive radius (Figure |22] | on dissipation, we would expect that 
systems with smaller at fixed stellar mass (an indicator of 
greater dissipation) should exhibit somewhat stronger stellar 
population gradients. Figure |32] shows the predicted gradi- 
ent strength versus effective radius from simulations and ob- 
served systems. We focus on the metallicity gradient: the 



30 



Hopkins et al. 



color and age gradients evolve strongly with time, introducing 
much larger scatter and possibly even reversing the predicted 
trends here. There is a significant trend predicted, and the 
observations show some tentative evidence for such behavior, 
but more data are needed to determine this robustly (because 
there is a mean offset in the range of metallicity gradients at 
different masses, it is important to consider narrow ranges in 
observed stellar mass). To the extent that the strength of the 
metallicity gradients is an indicator of the degree of dissipa- 
tion in the elliptical formation event, this prediction poten- 
tially provides a completely non-parametric means to search 
for non-homology (i.e. differential degrees of dissipation as a 
function of radius at fixed mass) along the fundamental plane. 

9. DISCUSSION 

We have studied the formation and properties of extra light 
and cusps in elliptical galaxies using a large library of both 
numerical simulations and local observed cuspy ellipticals. 
We demonst rate the relation of these ob served cusps to the 
prediction of iMihos & HernquistI (Il994ah with our improved 
numerical models: namely, that tidal torques in major merg- 
ers of gas-rich disks channel gas into the central regions of the 
galaxy, where it forms a dense central starburst. The starburst 
leaves a central light excess with a high phase space density, 
making the remnant more compact and reconciling the den- 
sities of disks and elliptical galaxies. Improvements in spa- 
tial resolution (the "extra light" in previous numerical studies 
has generally been entirely unresolved), numerical treatments 
(see Appendix iBl). and models for star formation allow us to 
make detailed comparisons and for the first time demonstrate 
good agreement between the extra light and cusps in simu- 
lated gas-rich merger remnants and observed cusp ellipticals. 

We argue that stars in these cuspy ellipticals (specifically, 
we intend this to represent those ellipticals formed immedi- 
ately in gas-rich major mergers; not those which may have had 
subsequent gas-poor re-mergers) should be separated into at 
least two distinct populations. First, stars which are formed in 
the disks (or otherwise in extended distributions in progenitor 
galaxies) before the final merger and coalescence of the two 
galaxies. The final merger scatters these stellar orbits and they 
undergo violent relaxation. They dominate the light, even 
in highly gas-rich merger remnants, outside of 0.5 - 1 kpc, 
and form a Sersic-law profile owing to their partial violent 
relaxation. Second, the starburst or dissipational population, 
formed in the central gas concentration in the final merger. 
This component is compact, and dominates the light inside a 
small radius < 0.5 - 1 kpc. These stars do not undergo signifi- 
cant violent relaxation, but form in a nearly fixed background 
potential set by the dissipationless component of the galaxy. 
The size of the dissipational component is set primarily by 
the radius at which it becomes self-gravitating; the gas is then 
generally stable against further collapse (even with cooling) 
and rapidly forms stars. 

9.1. Comparing Simulations and Observations: Empirical 
Decomposition of Light Profiles 

Observations and simulations have, in recent years, been 
driven to convergence on this multi-component description 

There is also a third component present in simulations but not prominent 
in light profile fitting: gas moved to large radii temporarily either by feedback 
or tidal effects, which settles into the relaxed remnant and re-forms small ro- 
tationally supp orted components (embedded disks, kinematical ly decoupled 
cores, etc.; e.g. IHemquist & Bameslll991klHopkins et al.l2008ct) . 



of elliptical light profiles, at least within the "cusp" pop- 
ulatio n (Kormendv 1999; Kormen dy et al.1 120081; ICdte et al.1 
l2006l l2007i; iFerrarese et al.ii2006) . In particular, the combi- 
nation of ground-based and high-resolution HST photometry 
has allowed observers, from a purely phenomenological per- 
spective, to recognize central extra light components beyond 
the inward extrapolation of an outer profile. The interpreta- 
tion proposed herein has been advanced as the explanation 
for these central profile shapes, supported in a purely empir- 
ical manner by gradients in e.g. stellar populations, isophotal 
shapes, and stellar kinematics. 

By applying observational methods to a large library of sim- 
ulations, we can directly compare simulations with different 
degrees of dissipation to observed light profiles, in order to 
determine what is required in a merger-driven scenario in or- 
der to reproduce the observations. We find that in all cases 
of ellipticals of moderate masses > O.IM*, we have simula- 
tions which provide excellent matches to the observed sys- 
tems (down to our resolution limits ^ 30-50pc), compara- 
ble to the typical point-to-point variance inherent in the sim- 
ulation surface brightness profiles (cr^ ^0.1). We can there- 
fore further use this ensemble of simulations to test different 
means of decomposing observed profiles in order to calibrate 
a method which reliably and robustly recovers the appropri- 
ate physical decomposition in those simulations. We demon- 
strate (Figures [T] & |2]i that an appropriately designed para- 
metric fit to the total light profile recovers, on average, the 
correct (known physically in the simulations) decomposition 
of the light - the mass fraction and effective radius of both the 
physical starburst component and, separately, the outer dissi- 
pationless component. 

The details of our methodology are discussed in § |4l es- 
sentially, we fit the surface brightness profiles to the sum of 
two Sersic components, an inner and outer profile. The quan- 
titative details of the decomposition are, of course, sensitive 
to the methodology and assumed functional forms; for this 
reason observations have often disagreed over the quantita- 
tive scales of extra light even for similar p rofiles (compare 
e.g. ICdte et al.ll2006t iKormendv etani2008l) . Fitting the en- 
tire light profile to a single Sersic or a core-Sersic law, for 
example, and comparing the central regions with the fit can 
yield physically meaningless values for the "extra light," even 
though the result is, in a global sense, a formally good fit to the 
profile. Physically motivated models and simulated light pro- 
files of the sort studied here are therefore necessary in order to 
make quantitative statements about extra light. Furthermore, 
a wide dynamic range is necessary; we find that the extra light 
is typically ^ 1 - 10% of the galaxy light and becomes impor- 
tant within ~ 500pc (~ Q.lRe). The extra light component, 
as it becomes larger, typically blends in more smoothly with 
the outer profile, and does not necessarily appear as a sharp 
departure from the outer light profile. Therefore, simultane- 
ous resolution of small scales ^ lOOpc where the extra light 
dominates, and large scales 20- 100 kpc where the shape 
of the outer, dissipationless component can be robustly de- 
termined, is a critical observational ingredient enabling these 
comparisons. 

We apply this decomposition to our simulations and to 
a number of observed samples: including the combination 
of detailed elliptic al surface brightness pro files over a wide 
dynami c range in iKormendy et al.l (l2008l) and iLauer et aP 
(l2007ah . and a large sample o f gas-rich merger remnants ob- 
served in iRothberg & J osepE (2004). The results from our 
fits are included in Tables [T] & |2] (see Figures l3lfT3] and Ap- 
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pendix |A|i. For the observed ellipticals, we consider only 
those which are observed to have central cusps in high- 
resolution HST imaging, as opposed to those with central 
cores or mass deficits, as the latter are commonly believed to 
form in gas-poor spheroid-spheroid remergers, which we do 
not model here. We instead undertake a comparison of those 
o bjects and simu l ations o f both gas-rich and gas-poor mergers 
in [Hopkins etan(l2008hh . 

Applying these methods to observations, then, provides a 
powerful new diagnostic of the formation history of galaxies. 
In the physical context of the models considered here, it is not 
only possible to estimate (from the existence of some central 
component) whether or not the formation of a given galaxy 
requires some dissipation, but to empirically, quantitatively 
estimate the degree of dissipation - i.e. "how gas-rich" was 
the history or progenitor of a given elliptical. We extend this 
to develop a number of new probes of galaxy merger history, 
and corresponding predictions and tests of the merger hypoth- 



9.2. Predictions and Observations 

(1) From either fitted extra light/dissipational component 
masses, or a direct comparison of simulations and observa- 
tions, we find that there is a clear trend of increasing dissipa- 
tion in the profiles of less massive ellipticals (Figures [T9ll20l l. 
The required initial gas fractions in mergers to form profiles 
that resemble the observed systems are a decreasing function 
of mass, and at all masses from ^ 10^- IO'^Mq span a 
range bracketed by the typical observed gas fractions of spiral 
galaxies at the same mass, at z = (bracketing the low end of 
the required gas fractions) and z ~ 2 - 3 (bracketing the high 
end). In terms of the mass fraction in the final dissipational 
or starburst component, the trend is even more clear (signifi- 
cance ^ 5 (t), and is reflected in the observed extra light com- 
ponents. The trend is given by Equation (O: 



(/st: 



1 + 
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We are able to, for the first time, observationally directly infer 
the degree of dissipation from the surface brightness profiles 
of observed ellipticals and use it to show that, in line with 
prediction s from e.g. their fundarnental plane correlations 
(lCarlber£|Q 986; GunnI [19871 iLakg IT989I: Ikormendvl IT989I: 
iHernquistet al.iil993") . dissipation is more important in lower 
mass elUpticals, reflecting the increasingly gas-rich nature of 
lower-mass progenitor spirals. 

(2) At a given mass, the degree of dissipation strongly in- 
fluences the size of the remnant. In both observations and 
simulations we demonstrate a tight correlation between effec- 
tive (half-light) radius at a given stellar mass and the inferred 
dissipational/extra light fraction (Figure |22] |. This owes to 
the compact nature of the central dissipational component - 
increasing the mass fraction in this component means that the 
half-light radius must be smaller This is also reflected as an 
increasing velocity dispersion with extra light fraction at fixed 
mass, but this effect is weaker (because of the role of dark 
matter in setting the total mass and velocity dispersion). 

This has important implications for the fundamental plane 
of elliptical galaxies - namely that dissipation is a key driver 
of systems on the fundamental plane (in terms of e.g. stel- 
lar mass, in optical bands there are obviously also substantial 
stellar population effects). Looking at ~ ellipticals along 



the sequence of increasing dissipational component and de- 
creasing effective radius, we directly show how this relates to 
homology-breaking in ellipticals (Figures I24ll25] l. Those with 
large effective radii and little dissipation at a given mass are 
well-described by a Sersic law (with our mean Sersic index 
n ^ 2.5-3) over the range ^ 0.01 -5/?^. Those with smaller 
effective radii and large dissipational components typically 
show an excess at small radii, reflecting the concentration of 
starburst light owing to dissipation. This grows with smaller 
radii and larger dissipational components, from prominence 
only at small radii <C 0.1 7?^ to ^ Q.5Re in extreme systems. 

We caution that there is not an obvious sense of the shape 
of this extra component - in some cases the whole profile 
could still be well-fitted by a single Sersic law (albeit with 
a higher index n^). In others, there is only a small excess 
in the ji-r plane, but it extends to large radii ~ Q.5Re and 
therefore contributes considerably to the mass. In yet others 
it is similar to more obvious extra light, and to that visible 
in our most extreme simulations, as a sharp departure from 
the outer Sersic law at small radii. In any case, however, the 
tight correlation between dissipational component strength 
and size of the elliptical can be directly seen to relate to a 
subtle non-homology in the central surface brightness of the 
galaxies, in both the observations and our simulations. 

(3) We find that the outer Sersic indices of cusp ellipti- 
cals are nearly constant as a function of stellar mass or any 
other properties (with median ~ 2.5-3 and scatter ^ 0.7; Fig- 
ures fT4lfT7] i. We emphasize that, given our two component 
decomposition, this outer Sersic index is only meaningful in 
the sense of reflecting those violently relaxed stellar popula- 
tions which were present in the progenitors before the final 
merger (fitting the entire profile to a single Sersic index can 
yield a very different result). This is because the dissipation- 
less component is relaxed under the influence of gravity, and 
is therefore roughly self-similar across scales. 

We note this constant Sersic index distribution is rele- 
vant for gas-r i ch mer ger remnants. As we demonstrate in 
iHopkins et al.l (l2008hh . subsequent gas-poor re-mergers will 
systematically raise the outer Sersic indices, and pseudob- 
ulges or other low-mass bulges may have systematically 
lower indices (and together these may drive a system- 
atic dependence of Sersic index on mass, owing to the 
cosmological dependence of formation history on mass). 
This has been seen in the Virgo cusp elliptical sample of 
iKormendv et alj ( l2008l) . and each of our observed samples 
independently confirms our prediction, and together they 
show no dependence of Sersic index over more than three 
orders of magnitude in stellar mass. Especially restricting to 
a sample of ellipticals which are observed in multiple bands 
and from multiple ground-based sources and considering 
only those results for which the fits (from different bands and 
instruments) do not disagree, we confirm both the predicted 
Sersic index distribution and lack of dependence on other 
galaxy properties. This agreement further suggests that cusp 
ellipticals are the direct remnants of gas-rich major mergers, 
without substantial subsequent re-merging. 

(4) Extra light fraction does not, in any predictive sense, 
correlate with the global kinematic properties (namely, ellip- 
ticity, isophotal shape an, and rotation (V /cr)*) of ellipticals in 
our simulations or in the observations. However, there is in a 
broad sense a change in these properties for large versus small 
extra light fractions (Figure l26b. At small extra Ught fractions 
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< 5-10%, remnants are more similar to t he known prop erties 
of dissipationless disk-disk mergers re.g. lBarnes|[T988l) . with 
a wide range in rotation properties, isophotal shapes, and 
ellipticities. At substantial extra light fractions, however, the 
systems are (on average) somewhat rounder, and are more 
uniformly disky and rapidly rotating. The dissipational com- 
ponent itself makes the remnant potential more spherical, and 
is closely related to kinematic subsystems formed from gas 
that survives the merger w hich contribute to the rotation and 
isopho t al sh apes (see e.g. Cox et al ."2006b'; Robertso n et al.| 
2006bl: iDeke l & Cox 2006; Naab et al. 2006a; Burkert et alj 
2003). A. distribution of dissipational/extra light fractions 
like that observed yields remnants with a similar distribution 
in each of these properties to observed cusp ellipticals, in 
striking contrast with dissipationless simulations that have 
historically not been able to match these properties. 

(5) Black hole mass appears to be roughly correlated with 
the mass in the extra light component, in both simulations and 
observations, but we show that this is largely a reflection of 
the tighter correlation between black hole mass and total host 
spheroid stellar or dynamical mass (Figure l27b. This is not 
surprising: if black hole mass most tightly follows the central 
potential or binding energy of the bulge, as argued from 
recent observations favoring a black hole fundamental plane 
([Hopkins et al. 2007a b; Aller & Richstone 2007), then this 
is largely set by the total mass of the system. The extra light 
mass represents a perturbation to the binding energy/velocity 
dispersion the system would have if dissipationless, and there 
is in simulations a weak residual trend along the expected 
lines, but it is sufficiently weak that the present observations 
are inconclusive. 

(6) The dissipation of gas into the central regions of ellip- 
ticals in a central starburst gives rise to significant gradients 
in the stellar populations (Figure |28]|. Since they are formed 
at the end of the merger, from gas channeled to the center 
of the galaxy after being enriched by star formation through- 
out the progenitor disks, the central stellar populations tend 
to be younger and more metal rich. The detailed gradient 
structure in our simulations corresponds well to what is typi- 
cally observed - with metallicity increasing and age decreas- 
ing relatively smoothly from to ^ 0.03 -0.1/?^. We do 
occasionally find anomalous behavior resembling various ob- 
served systems - e.g. non-monotonic trends in the stellar pop- 
ulations - but these are generally rare, as is also observed. 

At large radii, observations find generally weaker gradients 
- our simulations suggest that these are not driven by dissipa- 
tion in the merger, but rather reflect pre-existing gradients in 
the initial disks. To be conservative, we have not included 
them here, but experiments suggest that they are not com- 
pletely destroyed: it is well known, for example, that particles 
roughly prese rve their ordering in binding energy in mergers 
(lBarneslll992h . so given the trend in spirals towards less metal 
enrichment and older ages in their outer regions, we expect 
a smeared out version of these gradients to continue weakly 
into the outer regions of merger remnants, where the stars are 
primarily those violently relaxed from the initial disks. Over 
the range where observed gradients are typically fit, we quan- 
tify the gradient strengths in our simulations and find they 
are similar to those observed, as a function of system mass, 
age, and metallicity (Figure |29]|. Previous claims that merg- 
ers would only yield weak gradients in ellipticals ignored the 
importance of dissipation, which allows us to form systems 



with metallicity gradients comparable to the most extreme ob- 
served (d[Z/H]/dlogr - -0.8). 

In general, the gradients are not strongly dependent on 
mass, but do depend somewhat on the degree of dissipation 
(Figure |29]|. They also evolve significantly with stellar popu- 
lation age, as suggested by observations - not only do the stars 
become uniformly older, but as they age after the merger, the 
difference in mass to light ratio between young stars formed in 
the starburst and older stars becomes less significant, so in any 
light-weighted stellar population the younger starburst stars 
become (relatively) less prominent, and the gradient weakly 
washes out (Figurel28Tl. Of course, being driven by the central 
dissipational population, the gradient strengths are also corre- 
lated with central metallicity and with one another, albeit with 
large scatter (and for e.g. uniformly high metallicities or weak 
gradients, these correlations break down). All of these trends 
are in good agreement with present-day observed ellipticals. 

Furthermore, the central stellar metallicities of ellipticals 
(both their mean and distribution about that mean) as a func- 
tion of e.g. mass and velocity dispersion, are reproduced in 
our simulations, given only the assumption that the initial 
spirals have metallicities appropriate for the observed disk 
galaxy mass-metallicity relation - in fact, when the initial 
disk gas fractions are moderate (> 20%), we obtain this result 
from the self-enrichment tracked in our code (which domi- 
nates the final total central metallicity) even if the initial gas 
and stars have zero metallicity (Figure l29ll. Given that we si- 
multaneously reproduce the central metallicities, metallicity 
gradients, and surface brightness profiles of observed cusp el- 
lipticals, the claim that these systems have too much mass 
in m etals to be p roduced from the merger of local spirals 
dNaab & Ostrikeiil2(j07l) does not stand up, at least insofar as 
we are restricting to cuspy ellipticals. 

We also reproduce the typical observed color gradients 
of these ellipticals (Figures |28] & |30] |. We caution that 
while time evolution of even physical quantities such as 
the metallicity and age gradients is non-trivial, the color 
gradients are even more complex. Generally, though, systems 
have blue cores (are more blue towards their centers) for a 
short time after the merger, commonly taken as a signature 
of a recent merger or central starburst, owing to the young 
central starburst population. However, this population is also 
more metal rich, so the gradient rapidly weakens and after a 
short period, ^ Gyr, the metallicity difference becomes more 
important and the systems have red cores (are redder towards 
their centers). The gradients are fairly weak, a difference 
of '--^ 0.1 -0.2 mag between the color within 7?e/8 and Re, 
and the simulated distribution of color gradient strengths 
agrees well with that observed in samples of local ellipticals 
(especially if we restrict broadly to similar mean stellar ages). 

(7) Since increasing the amount of dissipation at a given to- 
tal mass tends both to yield smaller effective radii and stronger 
stellar population gradients, we predict that the two should be 
correlated (again, at fixed stellar mass - considering too wide 
a range in stellar mass will wash out this correlation, since 
the two scale differently in the mean with stellar mass; Fig- 
ure[32]l. The trend is not predicted to be dramatic, but should 
be observable, and comparison with recent observations does 
suggest that systems with smaller effective radii for their mass 
do have uniformly higher metallicity gradient strengths. Fu- 
ture observations should be able to test this more robustly, 
using metallicity gradients (we caution that the age gradients 
and, especially, color gradients are too age-sensitive to be con- 
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straining in this sense). 

To the extent that gas-rich progenitor disks (at least at low- 
moderate redshifts) are so because they have more extended 
star formation histories, this should also be reflected in the 
integrated stellar populations of the ellipticals. This is well 
known in an integrated sense - low mass ellipticals (which 
form from low-mass, gas-rich disks) tend to have more ex- 
tended sta r formation historie s with less g-eleme nt enrich- 
ment (e.g. lTragereTa l. 2000; Thomas et ani2005h . At fixed 
mass, however, if systems undergo mergers over roughly the 
same time period (which is generally true for systems of 
the same mass), then those with smaller effective radii (and 
lower Mdyn/M*) should have formed with more dissipation 
from more gas rich disks, which presumably had more ex- 
tended star formation histories. This should yield younger 
pre-merger stellar population ages and less a-enriched stellar 
populations, relative to progenitors of the same mass which 
had shorter star formation histories and exhausted their gas 
supply before the merger There is one important caveat - sys- 
tems of the same mass might also have had more gas-rich pro- 
genitors because they underwent mergers at very early times 
(making them ol der and more g-enric hed), but cosmologi- 
cal estimates (e.g.l Hopkins et al1l2008bh suggest that systems 
with such early mergers will usually have multiple subsequent 
mergers at later times, so they will grow significantly in mass 
and have their effective radii substantially modified by these 
additional processes (and may reflect the same trends as more 
typical systems by z ~ 0). In any case, this latter process is 
not dom inant at a given stellar mass. We study this in greater 
detail in lHopkins et alJ (l2008d) . 

9.3. Synopsis 

We have developed a paradigm in which to understand 
the structure of (cusp or extra light) ellipticals, in which 
there are fundamentally two stellar components: the relic 
of a a dissipational central starburst and a more extended 
viol ently relaxed component ( i ntroduced and discussed in 
e^. 'Miho s & Hernauisl)|1994ai: iKormendyl [1991 ICox et alJ 
Hooeb; K ormendv et al.l l2008l) . We have shown that these 
components can be separated with observations of sufficient 
quality, and can be used to constrain the formation histories of 
ellipticals and infer physically meaningful properties in a hi- 
erarchical formation scenario. This allows us to demonstrate 
that dissipation is critical in understanding the properties of 
ellipticals, including (but not limited to) the structure of their 
surface brightness profiles, their sizes, ellipticities, isophotal 
shapes and rotation, age, color, and metallicity gradients (and 
their evolution), and the gas content and properties of their 
progenitors. 

We have studied these properties and identified robust 
trends across a large library of simulations, in which we vary 
e.g. the galaxy masses, initial gas fractions, concentrations, 
halo masses, presence or absence of bulges, presence or ab- 
sence of black holes, feedback parameters from supernovae 
and stellar winds, orbital parameters and disk inclinations, 
and mass ratios of the merging galaxies. This range of param- 
eters allows us to identify the most important physics. Most 
of these choices, for example, affect the surface brightness 
profile, extra light mass and radius of the extra light, concen- 
tration and effective radius of the remnant, and even its ellip- 
ticity and isophotal shape only indirectly. Ultimately, what 
determines the structure of the remnant (insofar as the proper- 
ties we have considered) is, to first order, how much mass 
is in the dissipationless (violently relaxed) component ver- 



sus the dissipational/starburst component at the time of the 
final merger(s). Therefore, changing something like the or- 
bital parameters or initial galaxy structure can alter the rem- 
nant substantially, but predominantly only insofar as it af- 
fects the amount of gas which will be available at the time 
of the final coalescence of the galaxy nuclei (i.e. how much 
mass ends up in the starburst component, as oppo sed to be- 
ing violently relaxed in this final merger; see also ICox et all 
2006c). Moreover, merger-induced starbursts may not be the 
only source of dissipation (for example, stellar mass loss may 
replenis h the gas supply and lea d to new dissipational bursts, 
see e.g. ICiotti & Ostrikeill2007l) . and the merger history and 
series of induced dissipational events may be more com- 
plex than a single idea lized major merger (see e.g. lKobavashil 
'2004; Naab etalj2007), but for our purposes, all dissipational 
star formation will appear similar when observed and have 
the same effects (we are essentially measuring the integrated 
amount of dissipation). 

We have demonstrated that this makes predictions for how 
fundamental plane scalings arise, which we study further in 
Hopkins et al. (2008a). Given these decompositions, we can 
attempt to observationally test whether sufficient dissipation, 
as a function of stellar mass, could have occurred in the in- 
ner regions of ellipticals to explain the discrepancies between 
their central densities and those of their proposed progenitor 
spirals. We make a wide range of new predictions for the 
distributions of these properties and how they scale with the 
degree of dissipation, and how they should scale with each 
other and various other observational proxies for this degree 
of dissipation (which we define herein). We have predicted 
and shown (given these proxies) that dissipation is indeed 
more important (contributing a larger mass fraction) in low- 
mass ellipticals, in line with expectations based on how gas 
fractions are known to scale with disk mass. Testing all of 
these with better observations should be possible in the near 
future, with well-defined samples of ellipticals and continued 
improvements in mapping e.g. the surface brightness profiles, 
stellar populations and their gradients, and structural proper- 
ties of ellipticals over a wide dynamic range. 

To the extent that their parameter correlations and theoreti- 
cal formation scenarios are similar, these conclusions should 
also generalize to "classical" bulges in disk-dominated galax- 
ies (as opposed to "pseudo-bulges" formed in secular pro- 
cesses, which we do not model herein, and which may or may 
not show similar structure). We have excluded them from our 
observational study in this paper because the presence of a 
large disk greatly increases the uncertainties in profile fitting 
(and makes it difficult to robustly identify multi-component 
structure in the bulge), but emphasize that re-analyzing the 
SO and early-type spiral galaxies in our observational sam- 
ples demonstrates that they are in all cases consistent with 
our conclusions. Indeed, studies of the central profiles of 
classical bulges in SO-Sbc galaxies, in those cases with suf- 
ficient resolution and dynamic range and without much ob- 
struction from the outer disk, have begun to see evidence for 
central two-component bulge structure with extra light com- 
ponents similar to those predicted here (consistent in their 
profile shapes, sizes, mass fractions, kinema tics, isophotal 
shapes, stellar population s, and colors; see e.g. lBalcells et al.l 
120071; iPeletier et alJlIOOTi) . We emphasize that these are dis- 
tinct from the nuclear star clusters seen in some bulges, dis- 
cussed in Appendix ICl 

It is important to extend this study to central core ellipti- 
cals, increasingly believed to be shaped by subsequent gas- 
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poor (sph eroid-spheroid) re-merg ers, and we consider these 
objects in iHopkins et alJ (l2008hh However, since cuspy el- 
lipticals dominate the ~ population, most of the mass den- 
sity in present-day ellipticals is contained in those cuspy ob- 
jects that we study herein. We also strongly emphasize that, if 
core populations are indeed the re-mergers of initially cuspy 
ellipticals, then the extra light or dissipational components 
are not destroyed in these mergers. Indeed, most simula- 
tions suggest that the dense, central components are suffi- 
ciently tightly bound that they feel relatively little disturbance 
in a dry merger. Although black hole scouring may scat- 
ter stars from the central region and create such a core (e.g. 
iMilosavljevic et ani2002h . these stars represent only a small 
fraction of the extra light population. In other words, scour- 
ing will flatten the central power-law like behavior of cusps of 
extra light, but not fundamentally remove the ^ 500 pc central 
light concentrations remaining from a dissipational starburst. 
Furthermore, the impact of scouring will be restricted to mass 
scales of order the black hole mass, which is only ~ 10"^ 
the stellar mass, much smaller than the typically masses of 
the extra light compo nents inferred through our analysis (e.g. 
iHopkins et al.ll2008hl) . Given this, it should be possible to 
generalize our modeling and constraints to much more com- 
plex merger histories than the idealized single major merger 
scenarios considered in our simulation study, provided we 
recognize that the dissipational and dissipationless compo- 
nents are really integral sums over the dissipational and dissi- 
pationless events in the formation history of a given elliptical. 
Of course, subsequent gas poor mergers may have a number 
of other effects on the structural properties of the galaxies, 
modifying many of the properties we have argued are initially 
set by or correlated with the degree of dissipation - there- 
fore we have excluded core ellipticals from our observational 
samples, and wish to emphasize the importance of doing so in 
subsequent observational comparisons. 

As we noted in Paper I, there is considerable room for 
progress in modeling the extra light component itself and 
structure of nuclear profiles at small radii (^ lOOpc), where 
observations are making rapid progress; however, modeling 
these radii in a meaningful sense requires not just improved 
numerical resolution but also the inclusion of new physics that 
will allow simulations to self-consistently form the structures 
that would be resolved (giant molecular clouds, star forming 
regions and star clusters, individual supernova remnants, the 
multiple gas-phases of the ISM and their exchange, etc.). It 



is an important and ambitious goal that the next generation of 
studies move beyond the sub-resolution prescriptions neces- 
sary when modeling larger scales, and include realistic ISM 

structure important in galactic nuclei.^ 

Here, and in previous work (e.g. IHopkins et al .' 2006a'' d. 
, 2008 b f ) we have developed a conceptual framework for un- 
derstanding the origin of starbursts, quasar activity, supermas- 
sive black holes, and elliptical galaxies. In this paradigm, 
the various objects and phenomena we have examined are 
connected in an evolutionary sequence that is ultimately trig- 
gered by mergers between gas-rich disk galaxies. Previously, 
from observed correlations between supermassive black holes 
and properties of their hosts, we have argued that ellip- 
tical galaxies must have originated through a process that 
blends together gas- and stellar-dynamics. This is moti- 
vated by simple physical considerations implying that su- 
permassiv e black holes are m ainly assembled by gas accre- 
tion (e .g. iLvnden-Belll 11969*) in order to satisfy the SoltanI 
( 119821) constraint relating black hole growth to quasar activ- 
ity, and the notion that elliptical galaxies were put together 
by violent relaxation (e.g. Lynden-Bell 1967) acting on stars. 
The results of our analys is applied to recent merger prod- 
ucts ( IHopkins et al.ll2008dl) . c uspy ellipticals (this pa per), and 
elliptical galaxies with cores (IHopkins et al.ll2008hh indicate 
that the same blend of gas- and stellar-dynamics is further es- 
sential for understanding the structural properties of the galax- 
ies that harbor supermassive black holes. Together, this pro- 
vides further support to the idea that elliptical galaxies and 
supermassive black holes originate via a common p hysical 
process (mergers of gas-rich galaxies) and, given the SoltanI 
(1982) constraint and the discovery of ULIRGs, quasars and 
starbursts, respectively, as well. 
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Cusps in Elliptical Galaxies 



TABLE 1 
Fits to Cusp Ellipticals 



Name 
(1) 


Ref. 

(2) 


Nphot 
(3) 


M, 
(4) 


My 

(5) 


a 
(6) 


Re 

(7) 


e 

(8) 


10004/0 
(9) 


(v/o-) 
(10) 


NGC 0596 


2 


4 


10.89 


-20.90 


164 


3.24 


0.07 


1.30 


0.67 


NGC 0636 


3 


3 


10.80 


-20.83 


156 


3.81 


0.13 


0.80 


1.04 


NGC 0821* 


2 


2 


11.36 


-21.71 


209 


3.63 


0.38 


2.50 


0.70 


NGC 1172 


2 


1 


10.46 


-20.13 


113 


4.37 


0.09 


— 


— 


NGC 1199* 


3 


3 


10.95 


-21.49 


207 


3.23 


0.23 


— 


0.49 


NGC 1400 


3 


2 


9.88 


-18.45 


250 


0.67 


0.13 


0.00 


— 


NGC 1426 


2 


3 


10.83 


-20.78 


153 


3.63 


0.31 


0.01 


— 


NGC 1427 


2 


2 


10.83 


-20.79 


170 


3.16 


0.29 


— 


— 


NGC 1439 


2 


3 


10.85 


-20.82 


159 


4.07 


0.15 


— 


0.32 


NGC 2314* 


3 


2 


11.34 


-21.95 


290 


3.37 


0.18 




1.17 


NGC 2434 


2 


1 


11.14 


-21.33 


229 


3.63 


0.07 


— 


— 


NGC 2534* 


3 


1 


10.43 


-20.25 


— 


22.67 


— 


— 


— 


NGC 2693* 


3 


2 


11.58 


-22.59 


279 


6.65 


— 


— 


— 


NGC 2768* 


3 


2 


11.01 


-21.56 


198 


6.78 


0.68 


— 


0.28 


NGC 2778 


2 


1 


9.67 


-18.75 


166 


1.62 


0.20 


3.50 


2.00 


NGC 2974 


2,3 


5 


11.00 


-21.09 


143 


4.07 


0.30 


0.50 


1.54 


NGC 3193 


3 


2 


10.51 


-20.30 


205 


1.81 


— 


0.30 


0.80 


NGC 3309* 


3 


3 


11.39 


-22.24 


266 


6.67 


0.15 


— 


— 


NGC 3311* 


3 


3 


11.91 


-23.14 


210 


12.15 


0.11 


— 


— 


NGC 3377 


1,2,3 


6 


10.42 


-20.07 


141 


2.24 


0.50 


0.50 


0.72 


NGC 3605 


2 


3 


10.16 


-19.61 


120 


1.58 


0.34 


-0.90 


0.74 


NGC 3610 


2,3 


6 


10.93 


-21.16 


163 


1.45 


0.52 


2.50 


1.10 


NGC 3818* 


3 


2 


10.24 


-19.69 


206 


1.61 


0.39 


2.30 


0.93 


NGC 3894* 


3 


1 


11.15 


-21.70 


— 


4.63 


— 


— 


— 


NGC 3904* 


3 


2 


10.93 


-21.07 


215 


3.05 


— 


0.00 


0.36 


NGC 3923* 


3 


3 


11.36 


-22.00 


216 


5.92 


— 


— 


— 


NGC 3962* 


3 


1 


10.85 


-21.04 


211 


10.31 


— 


— 


— 


NGC 4125 


3 


2 


11.39 


-22.53 


229 


4.97 


— 


0.95 


0.93 


NGC 4318 


1 


1 


9.83 


-16.26 


101 


0.59 


0.34 


0.40 


— 


NGC 4382* 


1,2,3 


4 


11.26 


-22.43 


196 


8.43 


0.19 


0.59 


0.33 


NGC 4387 


1,2,3 


5 


9.95 


-19.11 


84 


1.46 


0.43 


-1.00 


0.70 


NGC 4434 




1 


10.59 


-19.55 


118 


1.44 


0.08 


0.44 




NGC 4458 




1 


10.36 


-18.89 


85 


1.34 


0.12 


0.41 


0.32 


NGC 4459 




1 


10.56 


-20.88 


168 


3.30 


0.18 


0.22 


0.96 


NGC 4464 




1 


10.26 


-18.40 


120 


0.55 


0.30 


0.40 


— 


NGC 4467 




1 


9.71 


-16.87 


67 


0.39 


0.32 


0.60 


0.35 


NGC 4473 




1 


11.13 


-20.82 


192 


3.14 


0.39 


1.03 


0.28 


NGC 4476* 




3 


10.25 


-19.85 


41 


1.57 


0.28 


— 


0.60 


NGC 4478 


1,3 


3 


10.29 


-19.77 


149 


1.15 


0.19 


-0.80 


0.95 


NGC 4486a 




1 


10.36 


-18.84 


154 


0.50 


0.45 


4.00 


— 


NGC 4486b 




1 


9.90 


-17.70 


200 


0.19 


0.40 


0.80 


0.59 


NGC 4494 


2,3 


7 


11.24 


-21.50 


155 


3.72 


0.15 


0.30 


1.24 


NGC 4515 


1 


1 


10.03 


-18.52 


90 


0.87 


0.14 


2.00 


0.47 


NGC 4551 


1,2 


5 


10.02 


-19.08 


100 


1.25 


0.32 


-0.70 


0.55 


NGC 4564 


2,3 


7 


10.53 


-19.66 


153 


2.75 


0.45 


2.20 


1.05 


NGC 4621 


1,2,3 


7 


11.26 


-21.41 


237 


6.51 


0.34 


1.50 


0.81 


NGC 4660 


2 


2 


10.46 


-19.41 


191 


1.05 


0.40 


2.70 


1.04 


NGC 4697 


2,3 


4 


11.23 


-21.49 


165 


5.13 


0.41 


1.40 


0.78 


NGC 4742 


2,3 


4 


10.30 


-19.79 


93 


1.24 


— 


0.41 


1.62 


NGC 5018* 


3 


2 


11.36 


-22.29 


223 


4.17 


0.25 


— 


— 


NGC 5127* 


3 


1 


11.15 


-21.87 


— 


8.32 


— 


— 


— 


NGC 5444* 


3 


3 


11.20 


-21.91 


221 


5.46 


— 


— 


— 


NGC 5576 


3 


3 


10.81 


-20.82 


187 


3.52 


0.30 


-0.50 


0.22 


NGC 5638* 


3 


3 


10.65 


-20.60 


159 


1.99 


0.08 


0.20 


0.72 


NGC 5812 


3,4 


4 


10.97 


-21.12 


204 


2.63 


0.05 


0.00 


0.52 


NGC 5831 


3 


3 


10.68 


-20.42 


166 


2.56 


0.17 




0.19 


NGL. J 54 J 


3 


3 


10.38 


—19.46 


251 




n 1 c 


0.80 


0.91 


NGC 6482* 


3 


2 


11.57 


-22.76 


287 


4.51 


0.27 






NGC 6487* 


3 


3 


11.54 


-23.11 




12.91 








NGC 7562 


3 


2 


11.42 


-22.27 


243 


6.01 


0.29 






NGC 7626 


3 


3 


11.67 


-22.62 


234 


8.81 


0.13 


0.11 


0.12 


VCC 1199 


1 


1 


8.97 


-15.52 


69 


0.17 


0.28 


0.05 




VCC 1440 


1 


1 


9.13 


-16.75 


59 


0.70 


0.17 


0.60 




VCC 1627 


1 


1 


9.11 


-16.43 




0.30 


0.20 


0.40 




VCC 1871 


1 


1 


9.27 


-17.33 


51 


0.61 


0.14 


0.80 





ris (fit) 
(11) 



fij (sim) 
(12) 



/. (fit) 
(13) 



fsb (sim) 
(14) 
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0.080ii 
0.050 
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-0.009 
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0.090 
031+0-003 
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0.077 

034+0-008 
"•"■-'^-O.QIO 

007+O-002 

"•""'-0.002 
0.150 

0.075 
0.028 
0.228 
0.028 
0.148 
355+0-114 

n 22^+0.003 
"•^^^-0.195 

0.338 
0.073 

071+0-007 
"•""-0.018 

0.254 
0.040!0;025 
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0.094+ 



0.402 
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0.404 
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TABLE 1 — Continued 



Name 


Ref. 


Nphot 


M, 


Mv 


a 


Re 




(v/aT 


n, (fit) 


«s (sim) 


/. (fit) 


/* (sim) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


(11) 


(12) 


(13) 


(14) 



UGC 10638 


2 


1 


11.90 


-22.66 


14.79 






4.83 


ESO 462-15 


2 


2 


11.99 


-22.83 


289 7.24 0.27 








IC 2738 


2 


1 


11.64 


-22.21 


23.09 






6.33 


A0147-M1* 


2 


1 


11.86 


-22.60 


11.34 






3.86 


A0160-M1* 


2 


1 


12.19 


-23.18 


21.72 






1.54 


A0189-M1 


2 


1 


11.46 


-21.89 


9.33 






3.69 


A0261-M1 


2 


1 


12.06 


-22.95 


22.91 






6.85 


A0419-M1 


2 


1 


11.40 


-21.79 


8.13 






2.11 


A0912-M1 


2 


3 


11.66 


-22.24 


9.55 








A1308-M1* 


2 


1 


12.34 


-23.44 


26.66 






4.60 


A1836-Mr 


2 


2 


12.29 


-23.34 


22.05 






4 23+0- 13 


A1983-M1* 


2 


1 


11.72 


-22.35 


7.41 






3.79 



2.53: 

3.13: 

3.09 



-0.61 
1+0.15 



-0.18 
o ^4+0. IS 

2.64+0-57 

2 ^0+0.02 
-0 27 

2 ';q+i:io 

3 Ao+O-TS 
00+0.45 

^■"°-0.05 



0.001 
141+0.140 

"■'^'-0.140 

0.024 
0.030 
0.116 
0.017 
0.071 

0.153 

065+0.011 

0.043 
n ni s+fl ooi 
^•Ul'*-o.ooi 

0.033 



0.' 
0.099 



+0.093 

0.090|ii 

o.o78!!j:!J« 

0.090 



0.090: 

0.098 



-0.043 
+0.096 



-0.061 
105+*"* 

077+'' "'^ 

0.073i|:«i 

0.105+01?? 



Compiled and fitted parameters for the confirmed cusp ellipticals in our observed samples. Columns show: (1) Object name. (2) Source for surface brightness 
profiles, where 1 MKormendv etaH <2003) . 2 4Lauer et alj <2007al) . 3 ^ jBender et al] fMSl) . 4 = 4Rothberg & JosepHj2004l) . (3) Total number of different surface 
brightness profiles in our combined samples for the given object. (4) Stellar mass [logMt/M0]. (5) V-band absolute magnitude. (6) Velocity dispersion 
[km/s]. (7) Effective (half-light) radius of the total light profile [kpc]. (8) Ellipticity. (9) Boxy/diskyness. (10) Rotation. (11) Outer Sersic index ns of the 
two-component best-fit profile. Where multiple profiles are available for the same object, we show the median and minimum/maximum range of fitted n, 
values. (12) Range of outer Sersic indices fit in the same manner to the best-fit simulations, at ? fa 1 — 3Gyr after the merger when the system has relaxed. 
(13) Fraction of Ught in the inner or "extra light" component of the fits. Where multiple profiles are available for the same object, we show the median and 
minimum/maximum range of fitted values. (14) Fraction of light from stars produced in the central, merger-induced starburst in the best-fit simulations (it the 
approximate interquartile range allowed). This list includes all systems morphologically classified as ellipticals in lKormendv et al] j20O8l) (all are E0-E4), or 
(where not in that sample) Lauer et al. (2007a) (all are E or E/BCG). 

* Systems with ambiguous or uncertain cusp status. These are systems for which different sources disagree on their cusp/core status, or for which observations 
of the central regions are unavailable/ambiguous but for which some other evidence (e.g. stellar populations, gas/dust content, or kinematics) suggest a gas-rich 
merger origin. We include them here for completeness, but our conclusions are insensitive to their inclusion/exclusion, and they are not generally shown in our 
analysis. 
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TABLE 2 
Fits to Recent Merger Remnants 



Name 


Ret. 


Nphot 


M« 


Mk 


cr 


Re 


e 




(v/o-)* 


lis (fit) 
























NGC34 


4 


1 


11.09 


-24.61 


201 


0.84 


0.11 


1.36 


1.41 


3.49 


NGC 455 


4 


1 


11.10 


-24.64 


234 


3.33 


0.21 


-0.41 


1.08 


2.02 


NGC 828 


4 


1 


11.40 


-25.36 


- 


3.51 


0.41 


1.70 


- 


3.60 
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As Table[T] but for the recent merger remnant sample of lRothberg & Josephl 120041) . Note that we show the A'-band as opposed to V-band absolute magnitudes. 
Systems marked (*) are excluded from our comparison in this paper owing to contamination from a central AGN. Systems marked (f) should be regarded with 
caution, as unrelaxed or prominent disk/bar features make our fits unreliable. Some of the best-fit simulation parameters (Mi(sim) and /jj,(sim)) are slightly 
different from those in Paper I owing to an expanded set of simulations, but it makes no difference for our comparisons. 
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Fig. 33 . — As Figure[3] but in log-log space. Surface brightness profiles are shown for cuspy ellipticals in the Virgo cluster. Open circles show the observations, 
from lKormendv et~al . (2008). These are the highest-mass cusp or extra light ellipticals in Virg(i51(r^ 2M, ). Top: Observed V-band surface brightness profile with 
our two component best-fit model (solid, dashed, and dotted lines show the total, inner/extra light component, and outer/pre-starburst component). The best-fit 
outer Sersic index, extra light fraction, and variance about the fit are shown. Middle: Colored lines show the corresponding surface brightness profiles from 
the three simulations in our library which con'espond most closely to the observed system (shown outside the gravitational softening length, ^ 30 pc). Dashed 
line shows the profile of the starburst hght in the best-matching simulation. The range of outer Sersic indices in the simulations (i.e. across sightlines for these 
objects) and range of starburst mass fractions which match the observed profile are shown, with the variance of the observations about the best-fit simulatiorCSl 
Bottom: Observed disky/boxy-ness (04) and ellipticity profiles, with the median (solid) and 25 — 75% range (shaded) corresponding profile from the best-fitting 
simulations above. Note that these are not fitted for in any sense. Figures [34l41l show the other cusp ellipticals in the sample, ranked from most to least massive. 

APPENDIX 

A. FITS TO THE SAMPLE OF KORMENDY ET. AL. 2008 



In Figure s[33]|4T] we reproduce Figures [3]fT2l but with profiles shown in log-log projection as opposed to r'/"* projection. 
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Fig. 36. — Lower-mass cusp ellipticals (~ 0.2-0.3M»). Our simulations reproduce the observed outer profiles and kinematic properties of such galaxies, but 
do not resolve the stellar cluster nuclei at small radii. The extra light recovered by our two-component fits therefore can be misleading at such low mass. (As 
Figure|2] but in log-log space.) 
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Fig. 37. — Additional low-mass (~ 0.2M«) cusp ellipticals. Our fits perform better in this case. (As Figure[8] but in log-log space.) 
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Fig. 38. — Additional low-mass (~ 0.1— 0.2M, ) cusp ellipticals, but in this case without prominent stellar clusters in their nuclei. In this case our parameterized 
fitting is not misled and we recover similar starburst fractions to our simulations. (As Figure|9] but in log-log space.) 
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Fig. 39. — Very low-mass cusp ellipticals (~ 0.03 — 0.1 M, ). Our simulations provide less good matches at these luminosities, where dwarf galaxies dominate 
the spheroid populatio n (e llipticals at these masses are very rare). Robustly resolving the extra light in these very small systems probably requires < lOpc spatial 
resolution. (As Figure [Tol but in log-log space.) 
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Fig. 42. — Effects of resolution on a remnant from a simulation of a gas- 
lich merger. We show profiles of otherwise identical simulations with a large 
dissipational starburst, but different gravitational softening lengths e , . Dotted 
line is the profile (extrapolated to r — > 0) of a typical "cusp" elliptical (log- 
arithmic slope dln//dln(r) 0.5 as r — » 0). Finite resolution limits mean 
our profiles will always artificially flatten within some radius: triangles at 
top show where each simulation profile artificially flattens below this slope 
(dln//dln(r) = 0.5); this tends to happen at our softening length Ri e,. Arrows 
are shown for each simulation at 2ej: at 2 — 3e, the profiles are converged 
within ~ 0. 1 mag (good for our purposes in this paper). At > 5 e,, the profiles 
and kinematics are fully converged. 

B. ARE GAS-RICH MERGER REMNANTS "CUSPS"? 
RESOLUTION TESTS AS fi ^ 

In Paper I, we conduct resolution tests and demonstrate that 
the primary quantities of interest here, namely the extra light 
fraction and outer Sersic index, are reasonably well converged 
for the mass range of interest given our typical < lOOpc spa- 
tial resolution. In that paper, however, the observations had 
comparable (or poorer) spatial resolution to our simulations; 
the HST observations of the nuclear regions of ellipticals, on 
the other hand, resolve extremely small scales ~ 1 - lOpc, 
well below our typical simulations. It is therefore of some 
interest to examine the behavior at small radii in our simula- 
tions. 

Figure|42]shows an example of the (sightline-averaged) sur- 
face brightness profile of a simulated major merger remnant 
as a function of numerical resolution, reflected in the gravita- 
tional softening length Ej. Clearly, at some point around our 
resolution limits, the profiles artificially flatten (according to 
the softening) and become flat (essentially creating an artifi- 
cial "core" at the center).'"^ It is clear that, as we increase 
our resolution, the profiles continue to rise towards smaller 
and smaller radii. In detail, outside of ~ 2 - 3 £5, we find that 
the profiles are sufficiently well-converged for our purposes 
in this paper. However, more subtle, detailed features in the 
galaxies such as the boxy or disky-ness of the isophotal shapes 
04/0 (typically a ^ 1% effect in the deviation of the shapes 
from ellipses) and central kinematics, converge less rapidly. 
Nevertheless we find that outside of ~ 5 e,, there are no mea- 
surable resolution effects in any of these parameters. 

Where do our resolution limits flatten the profiles into false 

Note that this is different from the behavior seen in earlier generations 
of numerical simulations such as those in Mihos & Hernquist 1 1994a). In 
those simulations, time integration inaccuracies at the highest densities led 
the stellar cusps to artificially contract to the spatial resolution set by the 
gravitational softening length, yielding sharp, compact "spikes" in the sur- 
face brightness profile. This did not affect any other aspects of the evo- 
lution and, indeed, if these spikes are smoothed on the spatial scale set by 
the condition that th e starburst component be self-gravitating, the results of 
Mihos & Hernauisll fl994a ) a gree in detail with those presented here and in 
Hopkins etalj<2008dl 
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Fig. 43. — "Cuspy-ness" of our simulations at reasonably resolved radii. 
Top Left: Distribution of brightness profile slopes (7 = -dln//dln(r)) in the 
remnants of gas-rich simulations, at the smallest radii where resolution limits 
do not completely flatten the profiles (measured at 2.0 e, and 2.5 e,). Top 
Right: Same, but slopes are measured from the observed profiles at the HST 
resolution limits for ea ch object. Dotted li ne in both panels is a fit to the ob- 
served distribution (see lLauer et aU2007al) . Bottom Left: Profile slopes from 
the simulations, measured at fixed (fractional) radius (0.02/?^ and 0.03/?,.; 
we include only simulations where these radii are > Se,). Bottom Right: 
Observed profile slopes measured at the same fractional radii (0.02 and 
0.03 Rf). Dotted line in the lower panels is a fit to the observed distribution 
at O.OlRe. Down to our best resolution Hmits, the simulated gas-rich merger 
remnants show typical steep "cuspy" slopes similar to those in the observed 
cusp population. 

cores? We check this by simpl y adopting the common obser- 
vation al definition of a "core" dFaber et al.lll997t : iLauer et al] 
l2007ah . namely where the logarithmic derivative of the sur- 
face brightness profile flattens below a threshold / cx r"'/^ 
(-dln/(r)/dlnr < 0.5). This is shown in Figure|42l On av- 
erage, our systems only significantly flatten into false cores 
at r < 1 Ej. In other words, nearly all of our gas-rich merger 
simulations, including those pushing our spatial resolution to 
< 20 pc, have "cuspy" nuclear profiles (by the observational 
definition) all the way down to the gravitational softening 
length. 

This approaches (and in the best cases overlaps with) the 
radii where the observed slopes are classified as "cusps" or 
"cores." We therefore directly check whether the slopes of 
our simulated systems are similar to those observed. We al- 
ready know that many of our simulations match observed pro- 
files down to 50 pc, from Figures [3]fT3l here we use our 
highest-resolution simulations and test whether the slopes at 
the smallest radii we can reliably say anything about are rea- 
sonable or not. We consider the distribution in logarithmic 
slopes of our simulations at 2 e,, the smallest radii where they 
are not strongly flattened by resolution effects, a nd compare 
with the observed nuclear slope distribution from lLauer et alj 
( l2007al) at the observational resolution limits. In a more phys- 
ically motivated manner, we compare the slopes measured at 
fixed radii relative to the effective radius ^ 0.02- 0.03 in 
both our simulations and observations (comparison at fixed 
absolute radius yields similar results), only including simula- 
tions (and observations) where these radii are well-resolved 
(> 3 times the resolution limit). In either case, the agreement 
is good. This strongly suggests that, down to ^ lOpc where 
the physics of e.g. individual star-forming sites becomes im- 
portant, our gas-rich merger remnants genuinely have "cuspy" 
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Fig. 44. — Convergence properties of the quantities of inte rest in tliis paper, 
for an otherwise identical simulation (shown in Figure I42t as a function of 
gravitational softening length es. Filled diamonds show results from fitting 
the projected profile to our two-component decomposition, with a free inner 
Sersic index/shape parameter nj; open circles adopt a fixed n' = 1. EiTor 
bars show the ±lcr range owing to shot noise, statistical fit degeneracies, 
and variation between different realizations. The quantities of interest here 
are well converged for our typical resolution limits ~ 30-50pc. For all 
resolutions (< lOOpc) where the quantities are converged, we obtain the same 
average answers using a free n' or fixed = 1 ; but when resolution or seeing 
is relatively poor (> 30pc), a fixed n'^ is useful to minimize the noise (simply 
owing to the limited dynamic range in the fits). 

or "power-law" central profiles, and that the agreement seen in 
Figures [3lfT3]between simulated and observed profiles would 
continue down to such radii if we only had the numerical 
resolution. Similar conclusions were obtained in ICox et all 
(|2006c), who also demonstrated that these conclusions (given 
a fixed amount of gas at the time of the final merger) are insen- 
sitive to the details of the numerical algorithm and feedback 
prescriptions in the simulations. 

Because the inner shape of the extra light clearly changes 
(continuing to rise towards small r) in Figure |42] we exam- 
ine the convergence properties in a similar resolution test in 
Figure |44] We plot the fitted extra light fraction /extra, effec- 
tive radius of the extra light component /?extra, effective radius 
of the whole elliptical inferred from the fit Rg, and outer Ser- 
sic index in otherwise identical simulations as a function 
of €.5. We show the median values across a large number of 
sightlines, but the sightline-to-sightline distribution behaves 
as a whole in the same manner (see Paper I for resolution tests 
demonstrating the convergence of the distribution of /extra and 
«s across sightlines). For each simulation, we generate alter- 
native realizations by randomly re-scattering the stars accord- 
ing to the smoothing kernel, and show the error bars corre- 
sponding to the ±1 CT range in fits to these realizations (unlike 
e.g. the variation from sightline-to-sightline, much of which 
is real in that it reflects actual asymmetries in the galaxy, this 
variation is purely a resolution effect, and should vanish as 
Cs 0). The error bars therefore effectively include the for- 
mal statistical errors and fit degeneracies as well. We show 
results both for a fixed inner n , = 1 (the Sersic index of the ex- 
tra light component itself, as described in § |4|i, and free inner 
n' (i.e. both assuming a fixed extra light component shape or 
fitting for the shape). 

Two broad rules of thumb consistently emerge from our 
resolution tests. First, regardless of the exact details of our 
fit methodology, the quantities of interest here, especially 
/extra, and Rf are well-converged for spatial resolution be- 
low ^ lOOpc (~ 1" at the distance to Virgo, a factor 10 
larger than the HST diffraction limit). We demonstrate in Pa- 



per I that the distribution of these quantities across sightlines 
is also reasonably well-converged below this threshold. In 
terms of general criteria, global quantities (such as 7?^) con- 
verge quickly (at a resolution ^ a couple hundred pc), fol- 
lowed by integral quantities such as /extra and quantities re- 
lated to the "outer" (large-scale) profile (n,). The effective 
radius of the extra light (/^extra) is more demanding, because 
of course the fitted size will not be smaller than the resolution 
limits (generally we find this is converged up to a smooth- 
ing length < 1/2 of the converged or "true" extra light size). 
For this reason, we do not reproduce observed systems with 
^extra < lOOpc (see Figure ET] and Figures [39ll?n i. but this 
is only relevant for a few of the very lowest-mass ^ 0.01 L, 
ellipticals in our sample. Most demanding, of course, is the 
detailed shape of the extra light component itself (as r 0), 
which Figures |42] and |43] demonstrate we are only marginally 
resolving in our highest-resolution simulations. We reach al- 
most identical conclusions regarding convergence if we repeat 
this study by artificially degrading the seeing in the observed 
profiles. Note that the absolute values here are for ^ 0. 1 - 1 L* 
ellipticals - in much larger, more massive systems (especially 
those with flat central cores extending to ~ 100-300pc), the 
resolution limits can be much less restrictive. In any case our 
resolution studies and experiments with observed profiles sug- 
gests we are not significantly biased in our estimates of the 
most important quantity here, /extra (or /sb). 

Second, Figure |44] explicitly demonstrates that our results 
are not changed whether we (for convenience) adopt a fixed 
inner extra light component shape when we fit to our sim- 
ulations, or leave it as a free parameter. In the mean, the two 
recover the same answer (so long as our resolution is below 
the lOOpc threshold needed to resolve the structures of in- 
terest in the first place). Unsurprisingly, when the resolution 
is extremely good, the results are most robust when we allow 
a free n, - this allows the fit the freedom to deal with small 
scale features and e.g. the broad range in inner profile "cusp" 
slopes (Figurel43Tl. Fixing in such cases (or in e.g. observa- 
tions with ~ 1 - 5 pc resolution) can produce a higher rate of 
catastrophic failure owing to the presence of small-scale fea- 
tures that are unimportant for the overall profile. For example, 
a strict inner exponential (n' = 1) implies that the logarithmic 
slope of the surface brightness profile goes to zero as r ^ 0, 
whereas the small-scale (<C 30 pc) profiles of cusp galaxies 
have a wide range of non-zero logarithmic slopes. 

When the resolution is sufficient then, one can (and should) 
free «^ and fit for the shape of the inner component as well 
as its radius and mass fraction (although, again, fixing it in- 
troduces no bias, just a higher failure/confusion rate). How- 
ever, Figure |44] also shows the expected behavior as the res- 
olution is downgraded: up to the limits where properties are 
converged, fits with a free inner component shape recover the 
same answer on average, but the uncertainties (realization- 
to-realization noise and fit degeneracies) grow rapidly (ow- 
ing largely to the poor resolution of the converged extra light 
shape, which gives the fit too much freedom to trade off ex- 
tra light and outer components). The result is that, as noted 
in § m our results for this paper are entirely unchanged if we 
re-fit all our simulations with a free inner n', but the scatter 
in predicted quantities increases significantly. When the res- 
olution is poor, then, especially if the data sets of interest are 
limited (i.e. one does not have so large a number of objects 
that very large scatter is not a problem), fixing the inner shape 
is a convenient assumption that greatly reduces fitting degen- 
eracies and, for the choice of inner shape ~ 1, introduces 



Cusps in Elliptical Galaxies 



47 



no significant systematic bias. Fortunately, the behavior in 
Figure |44] and our other resolution tests is reasonably simple, 
and yields a useful rule of thumb: it is appropriate to free the 
inner component and fit for the shape of the extra light with 
resolution better than ^ 20-30pc (ideally < lOpc), or equiv- 
alently with resolution such that the extra light size /?extra is 
resolved with at least ^ 5 - 10 resolution elements. Again, we 
find the same is true for observations (the observations used in 
this paper - at least for the intermediate and massive galaxies 
- easily meet this criterion, with ^ 10-20pc corresponding 
to ^0.1 -0.3" at Virgo, a factor of a couple larger than the 
HST diffraction limit). 

Our dynamic range at large radii is, of course, not signifi- 
cantly limited in the simulations (rather we restrict to a com- 
parable range to that observed). For a detailed study of the 
dynamic range requ irements at large radii for fitting e.g. the 
outer H j, we refer to lKormendv et al.] (l2008l) . We note though 
their conclusions (which we also find in limited experiments 
changing our sampled dynamic range) that for the observa- 
tions and simulations here, the dynamic range is sufficient so 
as not to introduce bias or significantly larger uncertainties. 

C. THE DISTINCTION BETWEEN STELLAR NUCLEI AND 
EXTRA LIGHT 



Many spiral galaxies and dwarf spheroidals show a central 
excess in their light profiles associa ted w ith dense stel lar nu- 



clei ( 


nuclear star clusters) (Phillips et al.lll996; Carollo etal. 


19971 


1 9981 ICarolloll 1 999^ .Matthews et al.ll 1 999; ,B6ker et al. 


2003 


20041 iSeth et al.1 


20061). In their analysis of the ACS 


Virgo data, this is what 


Cote et al.1 (120061) andlFerrarese et all 



(|2006) identify as the central "excess" component (largely in 
their dwarf spheroidal sample). In rough terms, this aspect of 
these nuclei is superficially similar to what we identify as ex- 
tra light. However, as we discuss in §|4]and has been demon- 
strated with detailed HST observations of these objects (see 
ej. Carollo 1999; Matthews et al. 1999; Bokeret al. 200^ 
IWilcher et al.ll20 06). on closer examination it is immediately 
clear that these nuclei are very different and physically dis- 
tinct from the extra light or starburst components we identify 
in the observations and model in our simulations. 

Figure |45] compares the parameter correlations of extra 
light components and stellar nuclei; specifically their stellar 
masses, effective radii, and effective surface brightness or sur- 
face mass density ^ M /{2tiR^). The sequence of stellar nu- 
clei is clearly distinct from that of extra light. The typical ex- 
tra light component (effective radius ~ 100-500pc) is ~ 100 
times larger in spatial extent than a stellar nucleus (effective 
radius ^ 1 -5pc) at the same surface brightness. The slopes 
of the sequences in either projection shown in Figure |45] are 
nearly perpendicular: for starburst/extra light components, we 
find le oc M~*^-^ (less massive systems are more dense and com- 
pact); for the stellar nuclei, Ig oc M'l^-^ (less massive systems 
are less dense). If we include the limited subsa mple of stel- 
lar nu clei with velocity d i spersio ns measured in iGeha et al] 
(l2002l) and IWalcher et al.1 mO^, we find similar results in 
fundamental projections involving a. 

In fact, it is well established that the parameter correlations 
of stellar nucle i are similar to those of globular clusters, not 
ellipticals (e.g. ICaroUol [1991 iGeha et al.1 l2002t iBoker et al.1 
l2004t IWalcher et al.l I2005L and references therein), and we 
show this in Figure |45] The extra light components, on the 
other hand, form a relatively smooth extension of the se- 
quence obeyed by ellipticals (in terms of the elliptical total 



half-mass radius and mass; see e.g. lKormendvll985bh towards 
smaller radii and higher surface dens ities. This is a natura l 
prediction of our models: in § [T] and iHopkins et al.l (l2008ah 
we argue that ellipticals and bulges are (when analyzed as a 
single entity) driven along these correlations by the properties 
of their starburst/extra light, beginning from a location in pa- 
rameter space (occupied by the least dissipational ellipticals 
with large radii and low surface brightness) similar to their 
progenitor disks. In short, dissipational/starburst components 
represent a smoothly rising excess at ^ 0.5 - 1 kpc from an 
outer dissipationless component extending to ~ 10- 100 kpc 
in classical bulges and ellipticals, with smooth associated gra- 
dients; stellar nuclei represent a sharp excess (being a dis- 
tinct object with a fairly steep internal density profile; see e.g. 
IWalcher et al.ll2005l) at ~ 1-5 pc, with very distinctive prop- 
erties. To the extent that stellar nuclei exist in some ellipticals, 
they would sit "on top" of the extra light that dominates the 
profile within ^ 100-500pc, and exist entirely below the ra- 
dius regime we model or fit. 

As discussed in § H] this is borne out by a detailed com- 
pari son between ou r fits and those in e.g. Cote et al. (200(^ 
and iFerrarese et all (12006), who fit multi-component profiles 
to identify stellar nuclei in the very small-scale nuclear re- 
gions of Virgo galaxies. The "outer" profile in their images 
is based on the HST ACS profiles, which extent to outer radii 
^ 1 kpc (almost entirely dominated by the starburst compo- 
nent of the galaxy), and their identified stellar clusters typi- 
cally dominate the light profile at very small radii < 0.01 Rg. 
This is akin to separating our "inner" component itself into 
multiple sub-components - i.e. a starburst stellar component 
that blends relatively smoothly onto the outer, violently re- 
laxed stars and an innermost nuclear star cluster 

Moreover, the mass fractions of these components are very 
different. Figure |46] compares the distribution of stellar nu- 
cleus mass fractions (relative to their host bulges) versus 
the mass fractions of extra light/starbursts. The characteris- 
tic stellar nucleus has a mass f raction ~ () .001 - 0.003 Mpai 
(we plot the distribution fitted in ICote et akl f2006.). but other 
studies find very similar di stributions; see e.g. ICarollo et al.l 
ll997HMatthews et &\\\\99'^ . compared to ~ O.lMgai in extra 
light. Because the star clusters are typically identified in much 
less massive galaxies, the discrepancy is even larger 4-5 
orders of magnitude) if we consider absolute masses. No 
amount of hierarchical merging, then (which would roughly 
conserve mass fractions in these components), would move a 
stellar cluster into the regime of extra light components. 

A number of other properties reveal the sharp division be- 
tween these populations. The stellar populations in the extra 
light (see § |8]l tend to be of similar (albeit slightly younger) 
stellar age to their hosts (~ 3- lOGyr), and are highly metal 
enriched (typical central metallicities ~ 1.5-3 times solar) 
and moderately a-enhanced ([a/Fe] ^ 0.2-0.3). Stellar nu- 
clei tend to be extremely young (ages < lOOMyr) an d have 
somewhat sub-solar metallicities (IWalcher et al.ll20()6h . There 
are also kinematic differences, with characteristically less ro- 
tation in nuclear clusters (Geha et al. 2002). 

Altogether, there should be no risk of confusing stellar nu- 
clei and dissipational components with detailed observations. 
It is worth noting that at the very lowest elliptical masses 
(M, < 10*^ M0, Mb > -17) and very highest stellar nuclei 
masses, the parameter sequences in Figure |46] approach one 
another, albeit still with opposite slopes. The most likely ex- 
planation is that this is just a coincidental overlap of their 
structural scalings, and in any case, true ellipticals and clas- 
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Fig. 45. — Fundamental parameter correlations of extra light components (open red stars and violet circles; style as in Figure [T4l . compared to those of 
stellar nuclei (nuclear stellar clusters or "central massive objects"; blue symbols). Left: Effective surface mass density versus stellar mass (a nearly identical 
plot is obtained if we consider surfac e brightness versus luminosity). Observed stellar nuclei parameters are from Geha et al. 1 2002, triangles), Boker et al] 
12004 . inverted triangles), and. Walcher et al.i 12005 . diamonds). The two classes of systems not only separate strongly in mass, but trace nearly perpendicular 
correlations with opposite physical senses. For compaiison, we plot the corresponding points for the entire elHptic al profiles (fi l led symbols) an d for globular 
clusters (pink symbols; X 's and crosses are Milky Way and NGC 5128 (Cen A) globulars from HaiTis 1 1996) and Harris et alj <200a) ;IMartini & Ho (2004), 
respectively). Extra light components form a continuous extension of the spheroid/elliptical population ("classical" stellar bulges also lie along this correlation 
when plotted); unsurprising given that we argue they drive the effective radii and central surface brightness of the ellipticals. Stellar nuclei, on the other hand, 
appear to follow (at least roughly) an extension of the correlations for globular clusters. Right: Same, but showing surface density versus effective radius of each 
component/population. Although the sequences may approach each other at the very lowest elliptical masses (where classical bulges and ellipticals are rare, and 
where we may have misidentified some dwarf spheroidals as ellipticals), there is a sharp division between their structural properties, masses, radii, and parameter 
correlations (also, stellar population ages and metallicities; see text). The two are in general easy to separate and should not be confused. 

[36l show examples where our parameterized fitting may fall 
victim to this misclassification. In general however, even at 
the lowest classical bulge masses and highest stellar nuclei 
masses, the other striking differences in the systems (their 
stellar populations and kinematics) remain distinct: there is 
no continuity or "intermediate" class between the two popu- 
lations in our data or the literature. 

It is in principle possible that there would be correlations 
between extra light and stellar nuclei, in th at both might 
be formed by dissipational processes (see e.g. 'Milosavlievi^ 
12004: Seth et al. 2006). Moreover, at some level similar 
physics may be involved in determining e.g. the competition 
between star formation and gravitational collapse that deter- 
mines their size-mass relations. However, given their sep- 
aration in stellar populations, and given the prevalence of 
stellar nuclei in systems such as bulgeless disks and dwarf 
spheroidals that have manifestly not experienced major dis- 
sipational angular momentum loss, they probably are not di- 
rectly coupled. 
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Fig. 46. — Distribution of the mass/light fraction of the galaxy bulge in 
stella r nuclei (dashed; from the fit to the distribution observed in Cote et al. 
I200fl) versus our fitted extra fight components (solid). Stellar nuclei have 
typical mass fractions 2 — 3 X lO^^Mbui, whereas extra fight is typically 
~ O.lMbui, a two order-of-magnitude difference. The gap between the pop- 
ulations is even larger if we consider e.g. absolute masses, effective radii, 
masses relative to the whole galaxy (as opposed to the bulge), or our best-fit 
simulation /sj, (instead of the fitted /extra, which has somewhat more scatter). 
The systems are more than just distinct; it is not possible to assemble typical 
extra light components from any amount of hierarchical merging/aggregation 
of stellar nuclei. 

sical bulges become extremely rare at such low masses (Uke- 
wise, stellar nuclei at and above these masses also become 
rare). It is also possible that, at these masses, some of our 
smallest ellipticals are really misclassified dwarf spheroidals 
(which have significant stellar nuclei; our fitting procedures 
might then mistakenly call that nucleus the "extra light") or 
that we are in these couple of cases finding a real stellar nu- 
cleus in an elliptical and accidentally calling it the extra light 
(as opposed to recognizing the larger, less dense starburst 
component, on top of which such a stellar cluster would sit 
at the center of the surface brightness profile). Figures [7] & 



TABLE 3 
Extended Fit Results 



Name 
(1) 


Moiph. 

(2) 


Ref. 

(3) 


lie (extra) 
(4) 


(extra) 
(5) 


n'^ (extra) 
(6) 


^J^e (out) 
(V) 


Re (out) 

(8) 


(out) 
(9) 


/. (fit) 
(10) 


«j (sim) 
(11) 


(sim) 
(12) 


NGC 4621 


E4/E/E5 


1,2,3 


i7.38rj:i 


-0.21 


a Q1+4.50 
-2.91 


-2.20 


■■^ -1.48 


a 9^+1.51 
-^-0.75 


053+''-'"'5 


3 06+'-™ 

-0.18 


n 140+0.168 



As Table[T] but with a complete list of fit parameters (included as supplemental material in the on-line edition of the journal; here, we show an illustrative example of one table entry). Parameters are the same as 
in Table[T] with additional fitted parameters listed. Added columns include: (2) Morphology (taken from the sources given in column (3), in the same order). (3) Source for surface brightness profiles (as in Table[T), 
where 1 = (Komiendv et alj 120081) , 2 = 4Lauer et al.l j2007af) . 3 = IBender et alj 419881) . 4 = tRothberg & JosephI 120041) . (4)-(10) Parameters of our two-component fits: (4) Effective surface brightness (in mag arcsec"' ), 
i.e. /x(Rextra), for the fitted inner (starburst) component. (5) Effective radius of the extra light component Rextra [kpc]. (6) Sersic index of the extra light component n'^ (cases with n'^ = 1 and errors = are where n[ is 
held fixed as described in the text). (7) Effective surface brightness of the outer (violently relaxed envelope) component. (8) Effective radius of this component. (9) Sersic index of this component («j(fit) in Table[T). 
(10) Integrated mass/light fraction in the fitted "extra" (inner) component. (11) Range of outer Sersic indices (equivalent of n5(out)) fit in the same manner to the best-fit simulations, at / 1 — 3 Gyr after the merger 
when the system has relaxed. (12) Fraction of light from stars produced in the central, merger-induced starburst in the best-fit simulations (it the approximate interquartile range allowed). These values are medians 
for the fits to the available photometric profiles for each galaxy; they do not, together, represent the best fit to any particular individual measurements. 



